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Abstract. We consider the defocusing cubic non-linear Schrodinger equation on general closed (com- 
pact without boundary) Riemannian surfaces. The problem was shown to be locally well-posed in 
H S (M) for s > i in [5]. Global well-posedness for s 1 follows easily from conservation of energy 
and standard arguments. In this work, we extend the range of global well-posedness to s > |. This 
generalizes, without any loss in regularity, the results in 5 18 , where the same result is proved for 
the torus T 2 . The proof relies on the I-method of Colliander, Keel, StafBlani, Takaoka, and Tao, 
a semi-classical bilinear Strichartz estimate proved by the author in |22l . and spectral localization 
estimates for products of eigenfunctions, which is essential to develop multilinear spectral analysis on 
general compact manifolds. 
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1. Introduction 



We consider the defocusing cubic non-linear Schrodinger equation on a 2D-Riemannian manifold 
{M,g a p) given by: 



idtu + A g u = \u\ 2 u 
it(0, x) = uq(x) 



(1.1) 
(1.2) 



where t G M, x g M, and u(t,x) is a complex-valued function on R x Af. A g is the Laplace-Beltrami 
operator on M (negative operator) with respect to the metric g a p. This equation is an important model 
in several areas of physics, including laser optics, plasma physics, and Bose-Einstein condensates (see 
[33]). We will be most interested in the case when M is closed, i.e. compact without boundary. 



1.1. Local theory: The state of the art of the local existence question for (11.11) depends on the 
manifold M under consideration. When M = R 2 , Strichartz estimates allow to prove local well- 
posedness in fP(R 2 ) for s ^ (see [T3] or [TT][3S] for a survey). The well-posedness for s > is 
described as subcritical, i.e. the time of existence guaranteed by the local theory depends only on the 
size of the initial data in iP(R 2 ). In contrast, the local well-posedness for s = is critical in the sense 
that the time of existence depends on the profile of the data. This difference stems from the fact that 
s = is the critical regularity for which the H s norm is left invariant by the scaling symmetry enjoyed 
by the equation: 



u(t,x)^ ju( — ,-). (1.3) 

The range s > is the scaling subcritical range and it determines the range of subcritical local well- 
posedness as well. Moreover, the equation is known to be ill-posed in a certain sense for s < (the 
flow is not uniformly continuous in any neighborhood of 0) [15] . 

The first answer to the local existence question on a compact manifold was provided by Bourgain 
in the case of the torus [2]. In this paper, he derives periodic Strichartz estimates for the linear 
Schrodinger equation and uses them to prove subcritical local well-posedness for scaling-subcritical 
equations, including (II. ip on T 2 for s > 0. The loss of e derivatives (for arbitrary small e) in the 
critical Strichartz estimates forbids one from proving local well-posedness in L 2 by arguing as in the 
M. 2 case. This problem is still open. As on R 2 , the equation is known to be ill-posed for s < [T4] | 15 | . 
As a result, one concludes that on T 2 , as on R 2 , the scaling subcritical range s > determines the 
range of subcritical local-wellposedness. 
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The situation changes if one considers more complex geometries for the compact manifold M 2 . In fact, 
one expects that, in contrast to the case of the wave equation, the infinite speed of propagation of the 
linear Schrodinger equation would cause the geometry to play a role in the local theory. This is the case 
when M is the 2-sphere S 2 . Here, Burq, Gerard, and Tzvetkov proved that (|l.ip is subcritically locally 
well-posed for s > and C 3 — ill-posed (i.e. the flow is not C 3 at the origin) for s < i cf. [2J[7]. The 
proof of well-posedness was based on a sharp bilinear Strichartz estimate which in turn is based on 
bilinear eigenfunction estimates and the sharp localization of the eigenvalues of the Laplace-Beltrami 
operator on S 2 . 

For general compact manifolds without boundary, Burq, Gerard, and Tzvetkov established in [5] (see 
also [35 ) Strichartz estimates for the Schrodinger equation on closed manifolds with a loss of derivative. 
This allowed them to prove local well-posedness of (jl.ll) for s > i. It is expected, as in the case of 
M. 2 , T 2 , and S 2 , that for any compact manifold M there exists a critical regularity sq(M) ^ (s = 
being the scaling critical regularity) for which one has subcritical well-posedness for s > sq(M) and 
ill-posedness for s < sq(M). The exact relation of this regularity sq(M) to the geometry of the manifold 
is a very interesting thing still to understand. 

We end our discussion of the local theory by mentioning that in the case when M is a compact manifold 
with boundary, local existence has been proved for all s > | in |31l building on the work of Blair, Smith 
and Sogge [5] . We should also mention the work of Ivanovici [35] on Strichartz estimates for compact 
manifolds with strictly concave boundary and exterior domains. 

1.2. Global theory: enjoys the following two important conservation law^J] 

• Conservation of mass 

M[u](t) := I \u(t,x)\ 2 dx (1.4) 
Jm 

• Conservation of the Hamiltonian 

E[u]{t) := J \\V g u{t,x)\ 2 + ~KM)| 4 dz - E[u}(0) (1.5) 

These conservation laws and the subcritical nature of the well-posedness for s = 1 imply global well- 
posedness in the energy space. Standard arguments imply global wellposedness for all s ^ 1 (see 
[35], [J]). The first result to prove local well-posedness of (|l.ip below the energy norm was due to 
Bourgain in the case when M = K 2 ([3][3])- Roughly speaking, the idea was to split the solution into a 
low frequency part supported on frequencies ^ N (N is a parameter chosen at the end depending on an 
arbitrarily chosen time interval [0,T]) and a high frequency one supported on frequencies larger than 
N. One then evolves the low frequency part by the nonlinear flow of (|1.1|) (which is global in time since 
low frequencies have finite energy) and the high frequency part by the linear flow e ltA (which conserves 
the H s norm). While the sum of those two flows is far from being a solution to (|1.1[) . the difference 
between the real solution u and this sum can be shown to be smoother, in fact in H , which means 
that the solution cannot blowup in 

This "high-low" Fourier Truncation method of Bourgain was the impetus to a more powerful method 
of Colliandcr, Keel, Staffilani, Takaoka, and Tao based on the almost conservation of a modified energy 
functional(see [16j for instance). Here one considers a Fourier multiplier I which is the Identity for 
low frequencies and an Integration by an order of (1 — s)(enough to make an H s function in H 1 ) for 
high frequencies. While the energy of the modified solution Iu is not conserved in general, it is almost 
conserved and this is enough to close an iteration argument and prove that E[Iu] does not blow up, and 
hence 1 1 it (t ) 1 1 jy « , being controlled by the latter, does not blow up either. The success of this strategy, 



There is also a conservation law for the momentum, but we will not be using it. 
2 The situation is actually a bit more complicated. One has to do this analysis on small intervals for technical reasons, 
which requires iterating the above-mentioned procedure to conclude that the interval of existence [0, T] can be chosen 
arbitrarily. 
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commonly referred to as the "I-method" for obvious reasons, is based on multilinear analysis in Fourier 
space and known local theory. 

In the context of compact manifolds, global well-posedness below energy norm was first considered by 
Bourgain in [5] where he used the language of normal forms. The problem was also studied in |T8] 
where a language similar to that on R 2 was used. The result in these papers is global well-posedness 
of JO) in H S (T 2 ) for all s > §. The proof relies heavily on a Fourier space analysis, which is one of 
the many advantages of working on an abelian group like T 2 or K 2 . We should also mention the recent 
paper [T] which proves global well-posedness on sphere-like manifolds (Zoll manifolds) in H s for all 
s > y| as well as some recent developments on the local and global theory for energy critical problems 
on compact, exterior, or product domains ( [251 [Ml 1251 1251 121)1 127] ) . 

In this paper, we consider the problem of global well-posedness below the energy norm for posed 
on a general closed (compact without boundary) Riemannian 2-manifold M. The first thing one misses 
when moving to the setting of compact manifolds is the Fourier transform. This loss is only partially 
compensated by the spectral resolution of the Laplace-Beltrami operator. Since the / operator is now 
given as a spectral multiplier rather than a Fourier multiplier as in the case of the torus T 2 (or on R 2 ), 
there are considerable difficulties in running the multilinear analysis in the former case. Firstly, in order 
to be able to perform a multilinear analysis with Littlewood-Paley pieces at fine scales, one has to know 
the spectral localization (which stands for the frequency support) of products of eigenfunctions. The 
resolution of this difficulty will be discussed below and leads to implications that we think are interesting 
in their own right. Secondly, after estimating individual products of eigenfunctions, one is faced with 
the problem of summing over these products without incurring additional losses in regularitjj^. This 
particular difficulty stems from the absence of a good analogue to Fourier inversion which allows the 
treatment of multilinear Fourier multipliers on K" or T n as convolution operators in the Fourier space. 
To resolve this latter difficulty, we use a Fourier series decomposition trick (cf. [13] [30]). The basic 
idea is to convert certain multilinear multipliers satisfying Coifman-Meyer type symbol estimates into 
tensored multipliers for which the problem can be easily solved because the sums decouple. This process 
is done using Fourier series expansions (see Section [5j. 

Another main difficulty encountered in the case of compact manifolds M in comparison to the torus, 
is the additional loss of derivatives in Strichartz estimates. The loss of \ of a derivative in linear (^in 
bilinear) Strichartz estimates not only leads to limiting local well-posedness to the range s > \ but 
also poses itself as a serious problem in establishing global well-posedness here as well. This should 
be compared to case of the torus where only e = 0+ derivatives are lost in the linear and bilinear 
estimates (see [2J,[S],[TH]) and to the slightly more suggestive case of the sphere S 2 where there is a 
necessary loss of \ of a derivative in the bilinear Strichartz estimate ([9]). Overcoming this difficulty 
is done by using the semi-classical bilinear Strichartz estimates proved by the author in |22j . In this 
latter paper, it is shown that at the semi-classical level, one can obtain bilinear improvement to linear 
Strichartz estimates similar (actually the same) as those refinements proved by Bourgain [3] on I 2 . 
While these short range refinements do not improve on the loss of half a derivative in the bilinear 
estimates on [0, 1] x M (at least without a smarter way of summing over the short range pieces of 
the interval [0, 1]), they do offer crucial improvements for bilinear estimates on the rescaled (inflated) 
manifold A M for A > 1. By moving the analysis to AM rather than M itself (for A growing large with 
the asymptotic parameter N), one is able to capture this gain and compensate (almost completely) 
the loss of derivatives on M itself in comparison to the case of the torus T 2 . 

After dealing with these new difficulties, we are able to run the I-method strategy and generalize 
without any loss in regularity the results proved in [5] [18] for T 2 : 

Theorem 1.3. Let M be any compact Riemannian 2-manifold without boundary, (jl.ip is globally well- 
posed in H s for all s > § . Moreover, the H s norm of the solution satisfies the following polynomial 
bound: 



'Any use of Weyl's law leads to a loss that cannot be tolerated. 
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\\u(t)\\ H s(M) <||«o||«. t^^ + . (1.6) 

1.4. Bilinear Strichartz estimates: We recall the bilinear Strichartz estimate on R 2 that was proved 
in [3] (and was first applied there in proving the above-mentioned "high-low" Fourier truncation 
method): Suppose that uq, vo € L 2 (M. 2 ) are frequency localized at dyadic scales Ni and N 2 respectively, 
i.e. suppu^ C : |£| ~ Ni} and supp-^ C {£ : |£| ~ ^2} with N 2 < iVi- Then 

fN 2 \ 1/2 

\\e lt u e lt vo||l 2 x (Rxr 2 ) < I — J ||mo||l 2 (r 2 )INHl 2 (r 2 )- (1-7) 

By a simple scaling and limiting argument, it is easy to see that fj 1 . T[) is equivalent to the short range 
inequality that is the same as CUT]) but with L^^Ni 1 ] x R 2 ) replacing Ll x {R x R 2 ). In [22], it 
is proved that this short range estimate holds true on any compact Riemannian manifold without 
boundary. The exact statement is as follows (From here on, we use the notation A to denote the 
Laplacc-Bcltrami operator A g ): 

Theorem 1.5. [22] Let M be a compact Riemannian 2-manifold without boundary. Suppose that 
uo,Wo £ L 2 (M) are spectrally localized on dyadic scales N\ and N 2 respectively, i.e. uo = l[Ni,2Ni){V~ A)uq 
and vo = lrjv 2 ,2JV 2 ) W~ ^) v o with N 2 ^ N\. Then 

{ N \ 1 / 2 

\\e ttA u e ztA v \\ L 2 {l ^_i_ ]xM) < (—) ||wo|Uj(M)ll«o||i»(AO- C 1 - 8 ) 
The same estimate holds if vq is supported at frequencies < N 2 . 

Since we won't prove this estimate in this paper, we remark that the numerology in (jl.8p can be 
understood (heuristically at least) by a simple back-of-the-envelope calculation. Thinking of e ltA uo as 
a "bump function" localized in frequency at scale N\ and initially (at t = 0) localized in space at scale 
■j^-. The Schrodinger evolution moves this bump function at a speed Ni thus expanding its support 

at this rate while keeping the L 2 norm conserved. Similarly, e ltA vo could be thought of as a "bump 
function" that is initially concentrated in space at scale ~ and moving (expanding) at speed N 2 . 
A simple schematic diagram allows to estimate the space-time overlap of the two expanding "bump 

N (d-l)/2 

functions" thus giving the estimate ^ 1/2 — for the Lf x ([0,N 1 ] x Mr) of the product. 

By splitting the time interval into pieces of length N± , one can use (|1.8[) to be get a bilinear Strichartz 
estimate over any time interval [0, T}. We will not be interested in the case of T = 1 (which is relevant 

for the local theory and for which the gain from the factor (^j^J is lost in the summatiorQ(see 

(|3.2jl )). but rather we will be concerned with estimates on small intervals [0, T] with T< 1 (actually 
iV-f 1 ^ T <C in most important cases). More precisely, our main use of Theorem 11.51 will be to 
overcome the difficulty of derivative loss in the bilinear Strichartz estimates on general closed manifolds. 
The idea to recapture the loss of derivatives is to "inflate" the manifold M and move the analysis to a 
dilated manifold M\ := AM with its rescaled metric g\. If A is large, one expects frequency-localized 
Strichartz estimates to become better behaved and exhibit a gain manifested by a negative power of A. 
This can be heuristically explained by the fact that, in the limit A — > 00, M\ becomes flat and hence 
the Strichartz estimates should resemble those on R 2 (where no derivative loss occurs). Theorem 11.51 
yields bilinear Stichartz estimates on [0,T] x M with (7W 2 ) 1/2 on the R.H.S. of ifO]) , By moving the 
analysis to the rescaled manifold M\ = AM (with A 1 chosen appropriately as dictated by a scaling 
argument in the I- method), this translates to a gain of A~ 1//2 in the bilinear Strichartz estimates on 



"^Unless one is able to find a smarter way to sum the contribution of each individual subinterval. 
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[0, 1] x M\ (cf. (|3.4|l ). a fact which allows to compensate completely for the extra loss of derivative on 
closed manifolds as compared to that on the torus T 2 . In short, at the scale relevant to the I-method, 
the Strichartz estimates deduced from Theorem 11.51 above are essentially equally "powerful" to those 
on the torus. 

1.6. Spectral localization of products of eigenfunctions: The I-method is based on a multilinear 
analysis of a multilinear spectral multiplier. Typically, in such analysis one splits all functions into 
Littlewood-Paley pieces and analyses each generic "product" by itself. One of the main reasons why 
frequency localization is such a useful tool in analyzing multilinear operators (ranging from products of 
functions fg to more general multilinear multipliers) on IR ci (or more generally on any abelian group) is 
the fact that if / is frequency localized in {|£| £ [N, 2N)} and g is frequency localized in {|£| £ [M, 2M)} 
with (say) N ^ M, then the product fg is frequency localized in the region { |£| £ [2N + 2M, N — 2M] . 
Ultimately, this is due to the fact that e 1 ^ 1,I e^ 2,1 = e l ^ 1+ ^ 2 '' x and the latter is also a character at 
frequency £i +£2- The same thing is true on the torus T d with ^1,^2 € Z d in which case the product 
of the two eigenfunctions e 1 ^ 1 ^ and e^ 2X is also an eigenfunction. 

Once we move to the case of general compact manifoldfQ, this sharp localization is not as straight- 
forward. Consequently, the following question becomes of importance: given two eigenfunctions f(x) 
and g(x) of the Laplace-Beltrami operator — A on a compact d— manifold M d with eigenvalues /z 2 and 
v 2 respectively (/x, v £ K+), what can be said about the spectral localization of the product func- 
tion f(x)g(x), in other words what can be said about 71^2 (/g) where tt„2 is the projection on the 
rj 2 — eigenspace of —A? For definiteness, we call this problem that of spectral localization of products 
of eigenfunctions. 

In the range when 77 ^> max(/i, v), an answer can be given using the parametrix representation of 
eigenfunctions and crude integration by parts. Essentially one obtains that | \ir n 2 (fg)\ \l 2 (m) ( r f)~ p 
for any p (see [5] and section 2]). However, this gives a very large range of localization of eigenvalues 
for the product fg, namely all eigenvalues rj 2 satisfying r] ^ Cmax(/i, v) for some large constant C. 

Studying the indicative case of the sphere suggests that the product should be localized at a much 
smaller range similar to that on the torus. On the sphere S d the eigenfunctions have a special form: 
they are given by spherical harmonics. Multiplying a spherical harmonic of degree k with another of 
degree of degree /, one can expand the product in terms of spherical harmonics of degree ^ k + 1. The 
eigenvalue associated with a spherical harmonic of degree k is k(k + d — 1) where d is the dimension 
of the sphere. This suggests that the spectrum of the product of two eigenfunctions / and g with 
eigenvalues /i 2 and v 2 respectively, and satisfying [i ^ v should be sharply concentrated around the 
range \/— A £ [fj, + u, fi — u] , in the sense that itrf-fg should decay when r\ ^ /i + v. 

Noticing that ^^(fg) is captured by the correlation integral J M h(x) f {x)g{x)dx where h is an eigen- 
function of —A with eigenvalue rj 2 , the problem is easily seen to be equivalent to obtaining decay 
estimates for such an integral. The most general result will be given in Theorem 14.21 and the remarks 
following it (where the more general problem of the spectral localization of a product of multiple 
eigenfunctions is studied). Theorem 14.21 provides an identity satisfied by the correlation integral (see 
equation (|4.3l0 . As a consequence of this identity one gets that if rj = \i + Kv with K > 1 (in other 
words if K := > 1) then the integral J M h(x)f(x)g(x)dx decays like K ~ J for any integer J (see 
Theorem 22] and remarks thereafter). The identity is proved using the Ricci commutation identities for 
covariant derivatives acting on tensor^ and an iteration argument. The proof is presented in section 
7. 

We leave the statement of Theorem l4.2l to section [U Here we cite a particular example of the estimates 
that can be deduced from identity (14.3[) . We note that the exact form of such estimates depends on the 



The results here apply to general compact manifolds M possibly with boundary and the eigenfunctions considered 
satisfy either Dirichlet or Neumann boundary conditions. 

®A11 needed notions from Riemannian geometry are revised in section [7| for completeness. 
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norms in which one would like to estimate the eigenfunctions or Littlewood-Paley pieces. For example, 
if one would like to estimate all eigenfunctions in L 2 (thus overruling the exponent distribution dictated 
by Holder's inequality), we have the following corollary that can be understood as a refinement of the 
bilinear eigenfunction (bilinear Sogge) estimates proved by Burq, Gerard, and Tzvetkov in [§I IIP]. 

Corollary 1.7. Let M d be a compact d— manifold possibly with boundary. Suppose that v 2 , A 2 ,/i 2 are 
eigenvalues of the operator — A with Dirichlet or Neumann boundary conditions satisfying v ^ A ^ /I. 
Ifu = X + K/j, + 2 Q/or some K > 1 (i.e. K = v ~ x - 2 > 1 ), then for any f,g£ L 2 (M) and any J € N: 



tt, (l [A , A+1] (V^A)/l [Al!A1+1] (V^A)g) <J -^P\\f\\L H M)\\g\\L^M) (1-9) 

where 

|y/ 2 4 /d = 2 

A(rf,M) := < /x 1/2 (log//)V 2 ifd= 3. (1.10) 

More generally, one has: 

J (v /z A)^l[A,A+i](v /T A)/l[ Al . Al +i](\/ Z A)5^ A ^ , / t ' ) ||fe||^(Af)ll/ll^(Af)l|g|l^(M)- 

(1.11) 

Interchanging the roles of A and v , the same estimates hold if v < A and ^ = A — if/i — 2 ?ot£/i > 1 
(i.e. if K >1). 

It is worth mentioning that the factor of A(d, fi) appears because we choose to estimate all eigenfunctions 
in L 2 and use a variant of the bilinear eigenfunction estimates of Burq, Gerard, and Tzvetkov|10j. In 
particular, the identity we establish in Theorem 14.21 does not exhibit any loss of derivatives, a fact 
which is often essential for performing an efficient analysis of the spectral localization of products of 
Littlewood-Paley pieces as well as products of eigenfunction clusters (see subsection 16 . 31 where we apply 
Theorem 14.21 to bound products of Littlewood-Paley pieces in a way that crucially does not lose any 
derivatives). 

The decay of K~ J mentioned above and guaranteed by Theorem l4.2l is enough to conduct the usual fine 
scale-course scale Littlewood-Paley interactions on a general manifold (even with boundary) similar to 
how it is done on R" and T n once one is able to overcome an additional difficulty mentioned previously 
and treated in section[S] In particular, it allows us to bound the notorious low-high frequency interaction 
in the I-method (see subsection 16. 3p 

The paper is organized as follows: in section [5] we fix the notation and explain our rescaling of the 
manifold M and its effects on eigenvalues, eigenfunctions, and linear Strichartz and X s ' b estimates 
which we recall as well. In section [l~4l we recall the bilinear Strichartz estimates proved in [52] and 
show that the estimates on the time interval [0,T] imply refined estimates on the rescaled manifold 
AM. We also prove a version of these estimates involving X s ' h spaces and with differential operators 
applied to the propagator e ltA uo- Such estimates will be needed after we apply the spectral localization 
machinery mentioned above and explained further in section|4l The proofs of the theorems in this latter 
section require recalling some ideas from Riemannian geometry and are delayed to the last section [7] 
in order not to distract the reader. In section [51 we deal with the problem of bounding multilinear 
spectral multipliers and formulate and prove a lemma pertaining to multipliers obeying Coifman- Meyer 
type estimates. In section we run the I-method strategy and prove Theorem II. 31 Finally, in section 
[3 we review the Ricci commutation identities and present the proofs of section U] 



the additive factor of 2 is used purely for technical reasons. 
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2. Preliminaries 



2.1. Notation: M denotes a closed (compact without boundary) C°° Riemannian 2-manifold (sur- 
face). In some sections (namely sections SI [5j and [7]) the analysis applies to general d— Riemannian 
manifolds possibly with boundary and M will also be used to denote such a manifold. We let g a p denote 
the Riemannian metric, g a/3 its inverse, V the induced (metric compatible torsion-free) connection, and 
A = V Q V Q the connection Laplacian. 

Throughout the paper, we use the notation X < Y to denote X ^ CY for some constant C > and 
A ~ B to denote A < B < A. All implicit constants are allowed to depend on M and on its dimension 
(in sections where we consider d— manifolds). Often times we will attach a subscript to < as in <a to 
denote the possible dependence of the implicit constant C on A. We also use the notation J+ when 
J € R to denote J + e for a fixed arbitrarily small positive number e. Similarly, J— refers to J - e. 
For functions / € Cg°(R x M) we denote by f(.,x) = J-tf(-,x) the Fourier transform in time of the 
function t H> f(t,x) given by f(r,x) = f R e~ ltT f(t, x)dt. We will often take the liberty of omitting 
the 2n factors in the definition of the inverse transform as these constants are inconsequential in the 
analysis. 

The Laplace-Beltrami operator on a compact manifold is a non-positive self adjoint operator with 
compact resolvenlQ. This gives an orthonormal basis of eigenfunctions {ek}^Li corresponding to a 
non-decreasing sequence of non-negative eigenvalues fik of —A. We will denote by Vf. the strictly 
increasing sequence of eigenvalues and by 7r„ fc the orthogonal projection on the Vk eigenspace. Often 
times (especially when we study the spectral localization of the operator \J — A) , it will be convenient 
to denote — n\ where £ R and 7r nfc = ir Uh . For any interval in / c [0,oo) we denote by Pj the 
orthogonal projection operator onto eigenvalues Vk with ^fvh S I, or equivalently Pj = A). In 

most cases, we will be interested in the case when / is a dyadic interval of the form [N, 2N) where 
N e 2 N in which case we denote P N = P[ N ,2N) for N = V with j = 1, 2, ... and Pi = P [0 ,2)- H S (M) is 
the natural Sobolev space associated with (Id— A) 1 / 2 with the norm: 



\\u\\h>{M) = ^{vk) s \U„ k u\\h(M) ) ~ E N 2s \\P N u\\ L 2 (M) J 

V k / \A'G2 H / 

We also define the space X s,b (R x M) as the completion of C^°(R x M) under the norm: 



\ V / - Vk) 2 \vk) S \\T^u{T)\\ L 2 {M) ) 

\ k T / 



l e ltAu \\H>>H°(RxM) 



where H^H^M. x M) = i? t b (R t ; H S S {M)) (similarly we use mixed Lebesgue spaces L q t L r x (R x M) = 
L?(R;L;(M)). 



®In section [4] and [7] we consider the Laplace-Beltrami operators on a compact d— manifold X with Dirichlet or 
Neumann boundary conditions. This also is an operator with compact resolvent on I? . 
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For every compact interval I C 1, we define the restriction space X s ' h (I x M) as the space of functions 
iioe/x M that admit extensions to R x M in X s ' b (R x M). X s > b (I x M) is equipped with the 



restriction norm: 



\\u\\x°,»(ixM) = inf {\\U\\x^»(RxM) with C/ = non I}. 

2.2. Rescaling Mi Let M be a C°° closed Riemannian surface. Any such surface can be thought of as 
being embedded in some ambient space R L . Using this embedding one can define a A-rescaled version of 
M for any A > 0, which we denote by AM =: M\ and is given by M\ = D\M where Da is the dilation 
by A in R L : x H> Ax. M\ inherits from R L (with its Euclidean metric) a metric ^0. When confusion 
might arise, we will distinguish tensors and operators on M\ with a A subscript or superscript, e.g. 
g*p,A\, etc... For a function / : M\ — > C, we will often denote / := D* x f the pull-back function 

given for y G M by f(y) = f(D x y) -■ /(Ay). It is easy to see that \\f\\ L p(M x ) = Ap||/|| L p(m) and 
Volume(M A ) = A 2 Volume(M). 

Since A\f(x) — j 7 Af(D^ 1 x) — p-(A/)(|)0 for any C 2 function / on M\, we get that the functions 
j e k( x ) form an orthonormal basis of L 2 (M\) with corresponding eigenvalues ^f-. As a result, the 
orthogonal spectral projection operator t^( v on M\ is related to 7r„ on M by the relation 



f(x) = * v foD?=(*j)fy. (2.1) 



2.3. Linear Strichartz estimates. We state some of the linear Strichartz estimates needed. These 
were obtained in [8] (see also [35]). We say that a pair of exponents (g, r) is 2d— Schrodinger admissible 
if2<g^oo,2^r<oo satisfy - + £ — ^. For any two admissible pairs (q, r), the following estimates 

hold for the linear propagator e ltA uo 



u o|U«Lj([0,l]xAf) ^q,r,M |K|| ff J (2-2) 



''/ii / ■-' / ; 1 1 i , l ; i ,, , A ; I Cl I 

' i (AT) 

(cf. 8 ). This estimate is derived from the following semi-classical Strichartz estimate: Suppose that 
uo = P[N : 2N)Uo, then 



l|e l u o|Uf£j([o,jv- 1 ]xM) < II u o||l2( M ). (2.3) 

Notice that this estimate has the same form as the corresponding estimate that holds on R 2 (where it 
holds on R x R 2 ). This theme will also show up when we treat bilinear Strichartz estimates in the next 
section. (12.21) follows by splitting the interval [0, 1] into N pieces of length each and applying (I2.3[) 
on each of them using the conservation of mass (a type of square function estimate is also needed, see 
i)- 



^From the differential topology point of view, the two manifolds M and M\ are the same. However, they are different 
from the point of view of Riemannian geometry. In fact, the rescaled manifold (M\, g^«) is isometric to the Riemannian 
manifold (M, X 2 g a p). The dilation map D\ : M — > M\ is a diffeomorphism from the manifold M to the manifold M\. It 
is crucial to note that the metric on M\ is not equal to {D\),g a p but rather satisfies the relationship D\ t g a p = -AjP^ o. 
The sample situation to keep in mind here is that of rescaling the interval [0, 1] to [0, A] . 

^This can either be seen from formula for the Laplacian in local coordinate given by r^—r- di \Z|3 A \9* 



_ |9 A I 

-At I — di \f\g\g 13 di since g x%3 = Yj-g lJ or by noting that g and \ 2 g induce the same connection V and hence 
A vlsl 
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Most of the analysis we will do is on the rescaled manifold M\. This requires us to know the dependence 
on A of the implicit constant in the linear (and bilinear) Strichartz estimates. This is possible by scaling 
and using the semi-classical estimate (|2.3[) rather than (|2.2[) . In fact, we will see here and in the next 
section that time 1 Strichartz estimates on M\ follow from time A- Strichartz estimates on M, which 
are obtained using (|2.3[) by splitting the time interval into pieces of relevant lengths. 

If ito G L 2 (M\) is spectrally localized on the interval [N, 2N), by which we mean that uq — PnUo, 
then the pull-back function uq e L 2 (M) is spectrally localized on the interval [AA, 2XN). As a result, 
we calculate: 



II HA X n \-|| it\ 2 A\ II \ 2 + -|| itA~ \\ 

\\ e u o||l«l-([o,i]xm a ) = X" \\e u \\ LqtL r i[0 ^ ]xMx) =Xi ''||e uo\\L q t L-([o,A I ]xM) 

since X 2 A\u (x) = (Au )(j )• If the length of the time interval js ^ then (|2.3p gives the bound 
(recall that | + | = 1): 



l|e l Au o||l«lj([o,i]xm a ) ^ M\ e * «olL?ij([o,^]xM) < M\uo\\l^m) - ||wo||l 2 (Mx)- 

If, on the other hand, ^ -rW, we split the interval [0, w] into y pieces of length and apply 
on each of them to get: 



fN\ 1/q fNX 1 ^ 
l|e l * a "o||l«lj([o,i]xm a ) ^ ( y J M\uq\\li(m) = f y J IKHl 2 (a/ a )- 

As a result, we get that: 



j (A M ^1 ||uo||l2(m) if A ^ A 

He -o|U ?W]x m,) < j (f) V« ||uo||i2(M) ifA<Ar (2-4) 

Notice that in the limit A — > oo, Ma becomes flat and the Strichartz estimates at a fixed frequency 
scale become the same as those satisfied on K 2 . 

We finally cite an additional Strichartz estimate that will be of use to us. Combining Bernstein's 
inequality for spectrally localized functions (which is true on compact manifolds, cf. Corollary 2.2 of 

8 

[8]) with the L\Lx semiclassical estimate in (|2.3[) one gets for u g L 2 {M) spectrally localized in the 
interval [TV, 2 A): 

Applying the same rescaling argument as the one used to get (|2.4[) we get: 

1 1 lt A ,, < J ^ 1/2 |ko||L2 (M ) if A^ A 

We u \\ L s J[0A]xMx) < j { Nf* N l/2 l[uollL2{M) iix ^ N ( 2 - 5 ) 

It is well known that any estimate of the form ||e ltA uo||y ^ II m o||l 2 where Y is a Banach space of 
space-time functions satisfying ||e z * 6 '/||>' < translates directly into an embedding of X 0,1 / 2+ into 
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Y (cf [50], lemma 2.9 of [35], or the next section for a bilinear version of this). As a result of this, we 
get the following consequences of (|2.4p and (|2.5|) : 



i j j ^ ) IMIx°.i/2+([ ,i]xj\/) ifA^iV 

y { — ) H w llx°- 1 /2+([o,l]xAJ-) II A iV 



iV 1 / 2 ||u|| x o,i/2 +([0 . 1]xM) ifA^iV 

|tt||£f t .([0,l]xJlf x ) ^ ) CATNl/8 1/2|| M , . ^-7J 

L IxJ iV H u llxo. 1 /2+([o,l]xM) ll A iV 



whenever u(t) = Pjvu(t). 



3. Bilinear Strichartz estimates 



In this section, we use the short-range/semi-classical bilinear Strichartz estimate (|1.8[) proved in [35] to 
derive bilinear estimates on [0, T] x M for any T > 0. These estimates over the interval [0, T] translate 
by rescaling to bilinear estimates on [0, 1] x M\. Since the spectral localization machinery of section 
Q] will eventually require us to bound products of the form P(D)e ltA uoQ(D)e ltA vo for differential 
operators P(D) and Q(D), we will have to obtain estimates for such products as well. Luckily, this 
follows directly from the parametrix representation of e ltA and the results in [32]. We also translate 
those bilinear Strichartz estimates into bilinear X s ' b estimates as was done in section |2~51 

The estimate (|1.8[) will serve as a building block for the estimate over the interval [0,1] on the rescaled 
manifold M\, In fact, by splitting the time interval [0,T] into ~ N\T pieces each of length ~ and 
using the conservation of L 2 norm one easily gets the following corollary: 

Corollary 3.1. Suppose the uq,vq 6 L (M) and N 2 ^ N± are dyadic scales. Then 

\\e ltA P Nl u e ltA P N2 v \\m[o,T]xM) < A(T,N 1 ,N 2 )\\ uo\\l 2 (m)\\vo\\l 2 (m) (3.1) 



where 



A(T,N 1 ,N 2 )< \ ($) ^fT^N 1 1 
[ (TN 2 ) 1/2 ifT^ TVf 1 



(3.2) 



The same estimate holds if vq is supported at frequencies < N 2 . 



Next we translate this estimate via scaling into a bilinear Strichartz esimate on [0, 1] x M\. Recall 
that if ek{x) is an L 2 normalized eigen-basis of the Laplacian of M with eigenvalues [ik, then \ek(j) 
is an L 2 (MA)-normalized eigen-basis of the Laplacian on M\ corresponding to eigenvalues ^f. As a 
corollary, we obtain the following: 

Corollary 3.2 (Time T estimate on M implies time 1 estimate on AM). Let Nx,N 2 G 2 Z and suppose 
u Q ,v e L 2 {M X ). If N 2 < Ni then: 



P Nl Uq e lt A A P N2 Vq 1 1 L 2 ( [0 ,1] x M x ) 



< 



< 



A(A- 2 ,AAri,AA^ 2 )||u ||L2 (MA) || t ; || L 2 (MA) 

1 /2 

(ff) IK||L2(Af A )IMlL2 ( M A) ifX^Nx 

(^) 1/2 I|uo||l2(m a )||wo||l2(a/ a ) ifX < Ni 



(3.3) 
(3.4) 
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The same estimate holds if vq is supported on frequencies < N 2 . 

Notice that in the "flat space" limit A — > oo, the bilinear estimate becomes the same as that on M 2 . 

Proof. The frequency localizations of uo and vo mean that ir\ (uq) — unless ^ G \Ni,2Ni) and 
tt\ ( Vo ) = unless e [N 2 ,2N 2 ). 

||e 4 * A ^ e JtA ^o||L 2 ([oa]xM A ) =A||e 4A2tA ^ e iA2tA ^o||L 2( [o,A-]xM A ) 

=A|| ]T e- l *( l ''= + ^)< VA2U o(^)7r A i/A2 «o(^)||L 2 ([oa- 2 ]xAf A ) 

where we denoted by irtjx 2 the orthogonal projection operator in L 2 (M\) onto the eigenspace corre- 
sponding to the eigenvalue vi/\ 2 . Define u : M — > C as u{x) = u(Xx). Using the fact that 

n v u(x) = ir* /x2 u(\x) 

we get: 



I iiA\ ziA\ 1 1 

\e x u e A vo||L2([ 0i i]xAf A ) 



=A|| £ e- <t ^ + ^7r 1 ^«o(fKt3b(|)|| i »( [ o,A-»]x^) 



and hence 



\\e UAx u e ltAx v \\ L 2 (%l]xMx) = \ 2 \\e itA uoe ltA vo\\mio,\-z]xM)- (3.5) 

Applying (|3.3[) we get: 

||e 4iA ^ e J * A ^o||L2 ([0il]xMA) <A 2 A(A- 2 ,A7V 1 ,AiV 2 )||uo|U=(M)||uo||L2(M) 

=A(A- 2 ,AAfi,A^ 2 )||u || L 2 (MA) ||w || L 2 (MA) . 



as claimed. 



□ 



Remark. The main gain provided by this lemma (at least for our purposes of applying the I-method) 
will be in the regime N 2 < A <C N\ in which case the second inequality in Q3.4)) will be crucial in order 
to get the full s > | well-posedness range. We should note that using the linear estimates alone is not 
sufficient even if one uses endpoint type estimates. In fact, in [31) . it is proved that if vq is spectrally 
localized at V^A e [N 2 ,2N 2 ) then: 

\\e ltA v \\ L 2 LrmN -i ]xM) < (logAr 2 ) 1/2 ||u |U2 (M) . 
Combining this with the trivial LfL 2 . and Holder one gets the following short-range estimate: 
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||e ltA -uoe rfA wo|| L 2( [0 jv-i ]xM) < (log A 2 ) 1/2 ||w ||l2( M ) ||u || L2(M ) (3.6) 
which would give at the time scale T: 

|, itA AtA„ || < j (logN2) 1 / 2 \\uo\\ L -(M)\\vo\\^(M) ^ T « N^ 1 

e u e Do y o,t]xJW ^ i ^1/2*^/2+ - f ^ > ^ (<*.7) 

I r 1 '-^' \\uo\\ L 2 (M) \\v \\mM) rfT>N 2 

which, in turn, would translate to the following estimate on My. 

H ua, lt A x I, ( [log(AAr 2 )] 1 / 2 ||u || i 2 (MA) ||uo|| L 2 (MA) ifA»A 2 

||e *u e a ^o||^([o,i]xm a ) < < w j/»t (3-8) 

^ ^t7^IN||l 2 (m a )IN||l 2 (m a ) it A < A 2 

When comparing (|3.4I) to (I3.8[) . one hrst notices a substantial improvement in the range A > Ni 

1/2 1/2 ... . 

versus (log AA 2 ) ). This is similar to the improvement provided by the bilinear refinement 



((f) 



to Strichartz estimate in [3] to linear Strichartz estimate on M. d . Another crucial improvement (espe- 
cially for our purposes of running the I- method) happens in the range A 2 <A< N%. In this range, 

1/2 

(|3.4p gives a bound of 1 < 1 whereas gives the large constant (log AA^) 1 ' 2 which is not 

enough to get the global well-posedness result. 



A couple of words about the proof of (|1.8[) seem to be in order. This will also allow us to justify 
a version of these bilinear estimates involving differential operators applied to the propagator e !t/ xj. 
The proof starts with the Burq, Gerard, Tzvetkov parametrix [5] of e ltA . Using this parametrix, one 
translates estimates like that in (11.8[) into bilinear oscillatory integral operator estimate of the form 
II^/'S'pSIIl 2 (KxR d ) for operators of the form: 



and 



T v f(t,x) = e^+^afrx&mdt (3.9) 



S lt g(t,x)= eW<t>*>Qb(t,x,Z)g(€)dZ (3.10) 



where u,fi > 0, a, b € Cg°(R x R d x R d ) and 4>,ip g C°° are real- valued phase functions satisfying 
a non-degeneracy condition and another crucial transversality condition (see [22] for details). This 
transversality condition is satisfied by the main terms of the parametrices considered at least when 
A 2 <C Ai. One then applies a bilinear oscillatory integral estimate (Theorem 1 of [22]) in order to 
obtain (j 1 . 8|) in the range Ai 3> A 2 . The case when A 2 ~ N\ follows from the linear Strichartz 
estimates. 

The parametrices of e N i u and e N 2 vq in [8] allow us to writd 1: 

e 7 ^ u (x) = T Nl u Q (t,x) + R Nl u (t,x) 



11 These are stated in Corollary 3.6 of 1221 
12 Strictly speaking this representation onl 
cover it by finitely many of such neighborhood, and hence we only need to prove the estimate on each one of them. 



12 Strictly speaking this representation only holds in an open neighborhood of xq £ M. Since M is compact, we can 
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and 



with 



e 1 ™ 2 A v Q (x) = S N2 v Q (t,x) + R N2 v (t,x) 



f Nl u (fx) = — W /_ j e iWl *(*' ,e 'Oo 1 (t,a!,€,JV 1 )fio(2ViOde (3-11) 



and 



S*Mt>*) = 7zh e iN ^ x ^a 2 (t,x^ 2 ,N 2 )v (N 2 ^)d^ 2 . (3.12) 

Here uq and #o are the respective microlocalizations of uq and vo in the considered coordinate patch 
(in particular ||«o||l 2 (r<<) < ||uo||l 2 (m) an d similarly for Vo)(ci. [5], [35]) and a%,a 2 £ Cq°(R x R d x R d ) 
are polynomials in and respectively. The remainder operators -Ra^ and i?Ar 2 are smoothing 
operators that satisfy: 



ll-RiVi^olUf H" ([-a.a]x M) N ± N \ \u \ \ L 2 (M) and | \Rn 2 V o | \hf H° {[-a, a] x M) <N A^ ^ I K I U 2 (M) 

(3.13) 

for any N. 

i - A 

If P(-D) is a differential operator on M of degree n, then P(D)e N i i*o has the following expression: 



P(£)e l — A u (x) = N?f' Ni u (t,x) + R' Nl u {t,x) 

where T' N and R' N are operators of the same form as and i?^ . In particular, T' N has an expression 

as in p. lip (just with different a) and R' N obeys the same estimates as in (|3.13p (by choosing N large 

• t 

enough) . Similar expressions for e N2 vq allow us to deduce that following from the exact same analysis 
used to prove (|1.8[) (see Corollary 3.6 of \'22\): 

Corollary 3.3. Suppose the Uq,vo G L 2 (M) are spectrally localized around Ni,N 2 G 2 Z respectively as 
in Corollary VS.'A Let P(D) and Q(D) be differential operators on M of orders n and m respectively: 



\\P(Dy tA u a Q(Dy tA v a \\ L 2 {[0 , T]xM) < 7V 1 "iV 2 m A(T,iV 1 ,iV 2 )|| W o||L 2 (Af)lko||L 2 (Af) (3.14) 
where A(T, Ni,N 2 ) is given in 



As before, we need to translate this "time T" estimate on M into a "time 1" estimate on M\. In order 
to simplify the scaling, we will only define define differential operators on M\ as rescalings of differential 
operators on M. Suppose P and Q are differential operators on M: for any / G C°°{M\) we define 
the operator P acting on / as P{D)f{x) = (D^ 1 )* o (P(D)f \ = (P(D)f \ (f ) where / G C°°(M) is 

defined as before by f(y) = /(Ay) for every y G M (In words, P(D)f is obtained by pulling / back 
to M to get /, applying P(D) to /, and finally pushing forward the resulting function to M\). With 
such conventions we have: 
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Corollary 3.4. Suppose the u n ,vn € L 2 (M\) are spectrally localized around N\,N% € 2 Z respectively 
as in Corollarv \3.1{ Let P{D) and Q(D) be differential operators on M of orders n and m respectively 
and define P(D) and Q(D) on C°°(M\) as indicated above. Then: 

\\P(D)e itA u Q(D)e itA v \\ LHlQA]xMx) < (AiV 1 )"(A7V 2 ) m A(A- 2 , A7V 1; AiV 2 )|| Uo || i2( M A )||«o||L 2 (M A )(3.15) 

< (XN^iXNzr i-fj \\u \\mM,)\\vo\\mM,) i/A < iVi(3.16) 

where A(T, Ni, N2) is given in Q3.2p . 

The proof is merely a rescaling of Corollary 13.151 performed as in the proof of Corollary 13.21 which we 
will not repeat. 

As for linear estimates, bilinear Strichartz estimates can be reformulated in terms of bilinear X S:b 
estimates. The reformulation of p. 31) and (13. 15)) is the following: 

Corollary 3.5. For any b > 1/2 and any f,g€L X°' b ([0, 1] x AM) spectrally localized in dyadic regions 
around Ni and N2 respectively (i.e. 1[jv 1 ,2JVi]('\/~ A)/ = / an d 1[jv 2 ,2JV 2 ] W~ A)<? — g), we have: 



ll/fl , IU? ia! ([0,i]xMx) ~ A ( A "^A^i, AA^ 2 )||/|Uo,i.([0,l]xAf A )l|g||x"-f([0,l]xM A )- (3-17) 

If P{D) and Q(D) are differential operators of orders n and m respectively defined as in Corollary \3.4\ 
then: 



\\P(D)fQ(D)g\\ L 2 J[as]xMx) < (AiV 1 )"(AiV 2 ) m A(A- 2 ,A^ 1 ,AiV 2 )||/|| x o, b([0 . 1]xMA) || 5 || x o, b([0il]xMA) . 

(3.18) 



Proof. We will only prove (|3.18l) as (13.17[) is merely a special case. Without loss of generality, it is 
enough to assume that /, g e Cg°([— 2,2] x M\). Let 



F(t) = e~ ltA f(t) and G{t) = e- ltA g(t) 



Then 

P(D)f(t) =P{D)e itA F{t) = f e itTl P{D)e itA F(T 1 )dn and 

Jr 

Q(D)g(t) =Q{D)e ltA G{t) = f e ltT2 Q(D)e ltA G(T 2 )dT 2 . 

JR 
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As a result, 

\\P(D)fQ(D)g\\ LHRxMx) 



f f e lt ^ +T ^P{D)e ltA F(T 1 )Q(D)e ttA G(T 2 )dT 1 dT 2 
P(D)e ltA F(T 1 )Q(D)e UA G(T 2 



L 2 ([-2,2]xAf A ) 

dT\dr 2 



L 2 ([-2,2]xAf A ) 

(XN 1 ) n (XN 2 ) m A(\-^,XN 1 ,XN 2 ) / / \\F(n)\\ LH M x )\\G(T 2 )\\ LH M x )dndT 2 
=(A7V 1 ) n (AAr 2 rA(A- 2 ,AAr 1 ,A7V 2 )||/|| x o, i >|| 5 || x o, 6 . 

□ 

We conclude this section with a statement of a standard trilincar X s ' b estimate that follows from 
p. 41) . In fact, it is well known (see [5] for example) that any bilinear Strichartz estimate of the form 
\\e ttA u e' ltA vo\\i J 2 < (A r 2 ) So ||Mo||L 2 ll u o||L 2 (where u and v are spectrally localized dyadicly around 
frequencies Ni ^ N 2 ) would imply a trilinear X s,b estimate of the form (|3.20[) . This is made precise in 
the following lemma borrowed from [§] (Proposition 2.5): 

Lemma 3.6. (Proposition 2.5 of Suppose that a bilinear estimate of the form 
||e ltA -uoe ltA wo||L 2 x ([o,i]xM) < (N 2 ) So \\u \\ L 2( M) \\v \\ L 2 {M) 

holds whenever u$ — lm lt 2N{\ (V — A)ito and vq = 1[at 2 . 2 7v 2 ] (\/—A)vo with Ni > N 2 . Then for any 
s > Sq there exists (b, b') £ R 2 satisfying 



< b' < - < b, b + b' < 1 



such that for any three functions iti, 1*2,1/3 £ X s,b (]R x M): 



(3.19) 



\uiu 2 u 3 \\ X3 , b > < \\ui\\ x ^i>\\u 2 \\x^b\\u 3 \\xs,b 



(3.20) 



The fact that IpTl]) holds on M A (since A(A~ 2 , XNi, XN 2 ) < (iY 2 ) 1/2 for A > 1) implies, with the same 
proof as in [5], that there exists (b, b') £ M 2 as in (13. 19)) such that the following estimate holds: 



\WlU 2 U 3 \\ X3 , b > ([Q1]xMx) < ||ui||x>.'»([0,l]xM x )||'"2|U=.l>([0,l]xM A )||'"3|U 3 . i '([0,l]xMx) ( 3 - 21 ) 

for any s > \ where the implicit constant is independent of A for A > 1. 

4. Spectral localization: Part I 

In this section, we deal with the problem of the spectral localization of products of eigenfunctions: 
given two eigenfunctions / and g of the Laplace-Beltrami operator on a compact d— manifold M, 
where is the product fg spectrally localized? In other words, what can be said about ir^fg) where 
■Kn is the projection on the fi— eigenspace of (—A). Of course, one can pose the same question for the 
product of any number of eigenfunctions. For our purposes, we will be most interested in the spectral 
concentration of the product of three eigenfunctions, but the same analysis carries on for any number of 
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eigenfunctions. The results for two eigenfunctions mentioned in the introduction can be easily obtained 
by setting the third to be the constant eigenfunction 1. We should also remark as well that the results 
in this section apply for any smooth compact d— dimensional Riemannian manifold M d including those 
with smootrFl boundary as long as one imposes either Dirichlet or Neumann boundary conditions on 
the eigenfunctions. 

As mentioned in the introduction, this problem is trivial in the case of the torus and the sphere because 
of special eigen-function bases in those two cases (characters e m ' x for T d and spherical harmonics for 
S d ). On a general compact manifold, the spectral localization of the product of two eigenfunctions 
e2 and e% with eigenvalues /12 and /J3 respectively on the /zi-eigenspace is detected via the inner 
products (ei, e2e3)L2(A/) where ei is an eigenfunction with eigenvalue fix. As a result, the above 
problem is reduced to that of estimating integrals of the form J M e\{x)e2{x)ez{x)dx with e\,e2,e% 
being eigenfunctions of —A with eigenvalues Hi,H2,H3- Of course, if one is interested in the spectral 
localization of a product of more that two eigenfunctions e2 and one needs to estimate integrals 
coming from (ei, {e.2 ■ ■ ■ efe)), namely j M e\{x)e2{x) . . .ek{x)dx. For the purpose of applying the I- 
method to the cubic nonlinear Schrddinger equation, we will be mainly interested in estimates for 
k = 4. 

The problem of identifying the spectral localization of the product of three eigenfunctions was encoun- 
tered in [5] where the following crude estimate was obtained: 

Lemma 4.1. Let ei, . . . be L 2 normalized eigenfunctions of the Laplacian with eigenvalues fi\, . . . /X4 
respectively. There exists C > such that, if for j — 2,3,4 we have C[ij ^ then for every p > 
there exists C v > such that: 



< <W- (4.1) 

A proof of this lemma can be found in [9] (lemma 2.6) and is based on the parametrix expression of 
the eigenfunctions and a simple non-stationary phase (integration by parts) argument. 

Unfortunately, the above estimate is way too crude for our purposes. To explain this point, we introduce 
the following notation. Since we will be interested in localizing — A g rather than — A g , it will 
be notationally convenient to denote, for i = 1,2,3,4, the eigenvalues fXi = nf where 7ij € [0, 00). 
We also assume without any loss of generality that n% ^ 713 ^ 77,4 and that 77,3 > (otherwise the 
answer is trivial). The suggestive cases of the torus T d and the sphere S d suggest that the integral 
J M ei{x)e2{x)e^{x)eA{x)dx should vanish (or at least present some sort of decay) if n\ > n2 + 113 + n± 
rather than n\ > C(ri2 + + 114) suggested by lemma [JT] The importance of such an improvement, 
in comparison to lemma 14. 1[ is most crucial when ri2 3> n% , in which case we are multiplying a high 
frequency eigenfunction 62 by two low frequency eigenfunctions e^ and e^. We expect that the resulting 
function, 626364, to be spectrally localized in the region n2 ± 0{n%) rather than the much larger region 
0(712) suggested by lemma |4~T1 

As a result, we are reduced to proving decay of the integral J M ei{e2e^ei)dx when n\ ^ 77,2 + -^77,3. We 
will be able to prove very fast (polynomial of any order) decay in terms of K of this integral. This will 
follow from the following theorem whose proof is postponed to the section [3 as it requires the revision 
of some ideas from Riemannian geometry. 

Theorem 4.2. Let M d be a compact d— dimensional C°° Riemannian manifold (possibly with bound- 
ary) and let e\, . . . , e^ be eigenfunctions of the Laplace- Beltrami operator (with Dirichlet or Neumann 
boundary conditions) corresponding to eigenvalues n\. ..n\ respectively. Denote by Aq: 



ia The exact regularity requirements will increase as a factor of the sharpness of the spectral localization one would 
like to prove. This is quantified by the number of iterations n in Theorem 14.21 which requires the eigenfunctions to be 
C 2 "(M) and hence it would be sufficient for the boundary (if it exists) to be C 2n . 



e\ {x)e2 (x)ea {x)e± [x)dx 



M 
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M 



e\(x) . . . e4,{x)dx. 



(4.2) 



Then for any n G N: 



A Q = 



-2) n A r 



(4.3) 



where A n is given by: 



A n = ei (-B„(e 2 ,e 3 ,e 4 ) + C„(e 2 , e 3 , e 4 )) cfa; (4.4) 

' M 



and B n (f,g,h) and C{f,g,h) are trilinear operators given by: 



B n (f,g,h) = O l+]+k=2n (V l f*Vg*\7 k h) (4.5) 

C n (f, .9, ft) = i+J - +fe <2(n-i) f-R» * V7 * V J ff * V k h) . (4.6) 

Here V 1 / denotes the i — th covariant derivative of f (which is a (i,0) tensor) and for any two tensors 
A and B , A * B denotes some contraction of A® B and (D(A * B) denotes a linear combination of 
contractions of A ® B . Here R n denotes a tensor obtained from the Riemann curvature tensor by 
contracting and differentiating it a bounded ( in terms of n) number of times. 

To see why this theorem provides the decay advertised above, we set n\ = n 2 + -KY13. This gives 
n\ - {n\ + n\ + nf) 2Kn 2 n 3 + (K 2 - 2)n§ ^ 2Kn 2 n 3 if K > y/2 and hence (|Q|l gives that 

I A) I ^ w ( n"l ) ■ Given the structure of A n in (14.41) . (14.51) and (14.61) each derivative falling on e 2 is 

J\ \ n 2 n 3 / 

accompanied by a factor of n 2 in the denominator and each derivative falling on e 3 or e 4 is accompanied 
by a factor of n 3 in the denominator, which makes the term ^"li essentially bounded at least from a 

heuristic point of vie w3. The proof of Thcorcm l4.2l is based on manipulations with covariant derivatives 
of the eigenfunction and is left to the last section [7] in order not to distract the reader. 

We should also mention that the fact that the operators B n and C n are multilinear allows one to use 
Theorem 14.21 which is a statement about single eigenfunctions e2,e 3 ,e4, to derive statements about 
eigenfunction clusters (see corollary below) or even Littlewood-Paley pieces of functions (see subsection 
16. 3p . A first instance of this is illustrated in the following proposition (cited in the introduction) which 
can be understood as a refinement of the bilinear Sogge estimates proved in j9| by Burq, Gerard, and 
Tzvetkov: 

Corollary 4.3. Suppose that A 2 ,/i 2 are eigenvalues of the operator — A g satisfying v A ^ /j. If 
v = A + K fi + 2 for some K > 1 (i.e. K = ^2^1 > 1 j ; then for any L 2 (M) and any J g N: 

7r„ (l [AjA+1] (V=A)/l^, M+1] (V=A) 5 ) <j ^fi\\f\\ L2m \\ g \\ L2(M) (4.7) 

14 The exact type of bound on ^"1 will depend on the context and the spaces involved. See for example Corollary 
14.31 An estimate that does not involve any loss of derivatives can be obtained if one uses Holder's inequality and estimates 
the eigenfunctions (or eigenfunction clusters) in L p spaces rather than L 2 . 



GWP OF CUBIC NLS 



19 



wher 



,1/2 



ifd = 2 

A(d,/x):=^M 1/a (logM) 1/2 ifd = 2, 

ifd^A 



(4.8) 



More generally, one has: 



J (\ /Z A)^l[A,A+i](\ /:i A)/l[ Al , Al +i](\/^A).9rfa; <j \\ h\\ L ^ M ) II /II l 2 (m) llfflU 2 (m) ■ 

(4.9) 

Interchanging the roles of X and v, the same estimates hold if v < X and v = X — K[i — 2 with K > 1 
(ic ifK:= >1). 



The proof of this corollary is also left to the appendix as it requires some ideas from the next section. 
It is essentially a consequence of Theorem 14.21 and the bilinear eigenfunction cluster estimates in [5] 
and [10]. It is worth mentioning that the spectral localization operator used here 1[a.a+i] W~ A) can 
be replaced by smoother versions like xW~~ A — A) with \ & <S(R). As before, this proposition is 
particularly useful when \i <C A. It says that that the product 1[a,a+i] W~ A)/ l^+i] (v 7 — A)g is 
spectrally localized (as measured by the L 2 norm of its projection onto various eigenspaces) in the 
region A + O(p) and starts to decay rapidly (faster than any polynomial power of K) as we move away 
from this region. 



5. A MULTILINEAR SPECTRAL MULTIPLIER LEMMA 

As is customary in previous applications of the I-method either on R d or T d (cf. [H],[T8]), one is faced 
with estimating k— linear multiplier forms: 



A(/i,.-.,/ fe )= / fh^i,^,...,^)h^i)-.-fk^k)d^i...^ 

J «i +...+& =0 ^ 

= J m (£i> • • • , Cfc-i, -6i - ••• - &-i) A(6) • ■ -fk{-£,i - ... - &-i)^i • ■ • Cfe-i 



One is then interested in proving L 2 or X s ' b type estimates for such forms, which, thanks to tools like 
Plancherel's theorem and Fourier inversion, can be done by proving weighted multilinear convolution 
estimates (in L 2 ) for (|5.1[) either on M. d or Z d . This is due to the fact that (|5.1[) is actually a weighted 
convolution in Fourier space. Such convolution estimates make sense if we replace M. d or Z rf by any 
other additive abelian group G and a systematic study of such estimates was done in |37j . 

Once we move to the realm of general compact Riemannian manifolds and away from the category 
of abelian groups, multilinear spectral multipliers cease to be expressible as multilinear convolution 
operators. In fact, the operators with which we will be concerned have the form: 



A(/i,...,/ fc ) 



E 

ni,...,n fc 



m(ni, ...,n k ) / % ni fi(x) . . 



Jk(x) dx 



(5.2) 



> M 



where fi(x) = ir m fi( x ) is the spectral expansion of / and 7r ni / is the projection on the nf—eigenspace. 
Note that in the case of the torus, one can use the Fourier expansion of / and Fourier inversion to 
write (|5.2p in the form (|5.ip with & G Z d . The estimates we will be interested in establishing take the 
form: 
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|A(/i,...,/ fc )| <« ||/i||y...||M|y 

where Y is some Banach space of functions like an L 2 Sobolev space or an X s,b space (in the latter 
case, one would be considering functions that depend on time and integrating in t as well as in x in 
the definition of A in ([521)). 

In the case of convolution expressions like (|5.1[) . one can use Fourier inversion to get such bounds. For 
example, in the case of the torus T d , one can bound (|5.2j) as follows: 



|A(/ 1) ...,/ fc )|<|M| L - / |/i(e l )|...|/*(C*)|de l ...de fc = IHU- / h(x)...f k (x)dx 

where fi(x) = |/(^)l e2 ■ One is then reduced to estimating the above space integral. This can 
be done by Holder's inequality for example (or any bilinear or quadrilinear estimate available). This 
would eventually give the desired estimate if one can show that 1 1 /| | y < ||/||y, which is trivial if Y is 
some L 2 based space like H s or X s ' h . 

One could try to follow the same strategy above in estimating (|5.2I) . Unfortunately, this leads directly 
to failure (except when k = 2) because taking absolute values |7r ni /i(a;)| forbids us to recover /j again 
partly because |7r ni /j(x)| is not an eigenfunction of the Laplacian anymore. 

The key to estimating expressions like (|5.2p for some bounded multipliers to is the following obser- 
vation. If to is given as a tensor product of functions depending only on one variable at a time, i.e. 
m(£i, . . . = mi(£i) . . . ?7ifc(£fc), one can postpone taking absolute values till after "inverting" the 
spectral decomposition. More precisely, 



k 



A(/i,...,/fc) = y~] I Y\{m(rii)ii n Ji{x))dx= I JT I y^mi(ni)it n Ji{x) \ dx = / fi(x) . . . f k {x)dx 
„. JM S , J Mi, \ „. Jm 



where fi(x) = Yli' m i{ n i)' K n i fi{x) is just a linear spectral multiplier applied to /. This would reduce 
the problem of estimating the multilinear multiplier into that of estimating the integral §Y\fidx in 
terms of 



If to is not a multiplier in tensor form, one can still estimate (|5.2p if to satisfies some smoothness and 
decay properties. The idea is to first split the rij sums into dyadic pieces rii ~ Ni and then use Fourier 
series to write fh{n\, . . . , 7ifc) = A{6)e % ( s ' Lni+ ~ + e * n <>) . This reduces to the case considered above 

since each e *( e i™i+ - + e fc™fc) j s obviously in tensor form (cf. [13] [30]). 

We won't formulate here the most general multilinear multiplier estimate in order not to distract the 
reader, but focus on a conditional lemma which is the case that will show up in the treatment of the 
I-method in the next section. For this, we assume that Y is Banach space of functions that satisfies 
ll/l|y ^5 ll/l|y whenever / = J2 n . e ieini Tr ni f is a frequency modulation of /. In this case, we will 
say that Y satisfies the "modulation stability property" . In particular, this is the case for X s,b spaces 
and more generally L 2 Sobolev spaces. Such spaces are the analogues of translation-invariant Banach 
spaces on M. d and T d . 

Lemma 5.1. Let A is a k-linear multiplier as in (|5.2[) associated with the multiplier m and let Y be a 
Banach space that satisfies the "modulation stability property" as above. Assume that fh satisfies the 



following symbol-type estimate^\: 
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(5.3) 



Suppose that the following estimate holds: 



o JM 



fi(t),...,f k {t)dxdt 



<S||/l||y...||/ fc ||y 



(5.4) 



whenever fi(t, .) = P^fift, .) are spectrally localized to frequency scales iVj. TTien £/iere exists a constant 
C ( depending only on implicit constants in (|5.4[) J smc/i £/ia£, 



A(/i («),..., <CB||/i||y...||/ fc | 



y ■ 



(5.5) 



Proof. Since n,; ~ iVj, one can write n,; = A^fij where rij € [0, 2). As a result, one can define a smooth 
function * € C c °°([-4,4]) such that: 



*(ni,...,n fe ) =m(iVini } ...,iV k n fc ) on [0,2) fc . 



(5.6) 



Thanks to (|5.3|) it is easy to see that has bounded derivatives of all orders (only finitely many orders 
are needed actually). Extending f as a periodic function to M fc allows us to express it as a Fourier 
series: 



# (fi l9 — , fife) = E A(0i,...,vix 



i(eifti + ...+0 fc n fc ) 



As a result of this, one can express A as follows: 



A(0 l5 ...,^)= V V A(0 1 ,...,0 fc )e 4(ei ^ + - +efc * ) / ■K n J i {t)...-K n J k {t)dx 



e a^,...,^) / f e 



,1/ 



E ^awU 



\n k ~N k 



= E A(6 u ... t 9 k ) [ ft 1 (t > x)..J e k "(t,x)d 



where /?' = X) n ~w e * ZNi7T nifi(t) for 1 < i < fc. As a result, 



15 



As is the case with usual multiplier theorems, one does not need to include derivatives of all orders for the theorem 
to be true. For this lemma, symbol estimates up to second derivatives are enough. 
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A(h,...J k )dt < Yl \MSu---,O k )\ I I f e 1 1 (t,x)...f e k "(t,x)dxdt 
J2 \A(9 1 ,...,6 k )\\\f e 1 i\\ Y ...\\f e k »\\ Y 



<dB J2 \M^---^k)\\\fi\\Y...\\fk\\Y^CB\\f 1 \\ Y ...\\f k \\ Y 

6»iGZ/4 

where we have used (|5.4j) in the second inequality, the stability of the the Y norm under spectral 
modulation in the third inequality, and finally the L 1 summability of A which comes from the fact that 
* is C 2 ([-4,4] fc ) (see for example [21]). □ 



6. The I-Method 

In this section we prove Theorem 11.31 by applying the I-method machinery. We first notice that if 
U(t,x) solves (|1.1[) on M over the interval [0,T], then the function: 



will solve the cubic nonlinear Schrodinger equation posed on the rescaled manifold M\, that is: 



(6.1) 



\u\ 2 u 



idtu + A\u = 

u(0,x)= u (x)eH s (M x ) 



(6.2) 
(6.3) 



over the interval [0, A 2 T0. We will be most interested in the range where < s < 1. Let N ^> 1 be a 
fixed large number to be specified later (depending only on T) and denote the eigenvalues of — by 
n\ < n 2 . < . . . where rn G K+ (Recall that nf = ^ where Vi are the eigenvalues of — A g ). We define 
the spectral multiplier: 



Iu = 22m(rii)-!r ni i 



(6.4) 



where m(rii) is given by: 



m{k) 



N\ 1 - s 



(f) 



if k < N 
if fc ^ 27V 



(6.5) 



and a smooth interpolent in between. For technical reasons, it will be preferable to specify m(fc) 
m (jj) where m : R — > [0, 1] € C°° is non-increasing and satisfies: 



m (t) 



1 if t ^ 1 

t-^ 1 - 5 ' if t > 2. 



(6.6) 



^We may assume without loss of generality that the initial data Uo is in C°°(M). One can remove this assumption 
after proving the polynomial bound on E using a standard limiting argument. 
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With this convention, Iu = J2n TO o(]v) 7r n i (u). Notice that / is the identity operator on frequencies 
^ N and is an integration operator for frequencies ^ 2N. As a result, I is smoothing of order 1 — s, 
in fact for any sq £ K: 

IMIh*o(m) < \\Iu\\ H s +(i->) < iV 1- *||«||H» . (6.7) 
In particular, for u £ H S (M\), Iu £ H l (M\) and we are justified to define the modified energy: 



E[u) = E[Iu] = [ \\\7 g Iu(t,x)\ 2 + -Uu^^x^dx. (6.8) 
Jm x 2 4 

Notice that E[u(t)] controls ||u(£)||^ s and hence the boundedness of E[u] over the interval [0, A 2 T] 
implies that ||u||jj s remains bounded (in particular it does not blowup). Our goal is then to show that 
for any T > the modified energy (j6.8[) remains bounded. 

Despite the fact that E[u] is not a conserved quantity in general (since Iu does not solve (|6.2j) in 
general), it is almost conserved: in the sense that its rate of change will be a negative power of TV. 
This will allow us to prove polynomial (in T) bounds for E[u\. 

The proof is in steps: First we prove that the equation satisfied by Iu is locally well-posed in X 1,b ([0, <5] x 
Mx) with 6 — |+ and S > 1. This will allow us to obtain the bound ^||/m || h i 1 (compare to 

the bound from (|6.7p obtained from the local well-posedness of (|6.2[1 ) that will be used in bounding the 
increment of the modified energy over the interval [0,(5]. The latter will be proved in Proposition 16.41 
using a multilinear analysis of dtE[u]{t). Finally, the process is repeated over sub-intervals of length 
5 partitioning [0,A 2 T] as long as the E[u](t) < 1. This latter condition will specify the number of 
iterations allowed and hence the dependence of N on T as well as the polynomial bounds on E[u(t)]. 

Remark. In contrast to previous applications of the I- method (cf |16j . [5].|18]). the increment on the 
modified energy will be bounded not only by a negative power of N but also by a negative power of 
the scaling parameter A (which will eventually be chosen to be a positive power of iV). This is the 
main reason why one is able to prove that global well-posedness holds on the full range s > 2/3 similar 
to that on the torus where better Strichartz estimates hold (but with worse dependence on A). The 
reason for this gain is the bilinear Strichartz estimates (|3.4p and the fact that A itself will be chosen in 
the end to be a positive power of N (namely A ~ N~^). 

6.1. Local well-posedness of the I-system: If u satisfies (|6.2p . then Iu satisfies: 



id t Iu + A g Iu= H\u\ 2 u) (6.9) 

lu(0,x)= Iu (x) £ H\M X ). (6.10) 

We shall need a local well-posedness result for the above I-system: 

Proposition 6.2. Suppose uq £ H S (M\) is such that \ \Iua\\jji < 1, then there exists < 5 ~ 1 such 
that (|6.9[) is locally well-posed on [0,<5] and \\Iu\\ x - Lt i + M ^ < 1. 

Proof. Iu satisfies the following integral equation on the interval [0, 5]: 



Iu(t) = e ltA Iu - f e ^- s)A I (\u(s)\ 2 u(s)) ds 
Jo 

and hence we have for some < b' < i < b with 6 + 6' < 1 : 
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II^IUlHIO.SJxMx) < \\ Iu o\\hHM x ) +S 1 (6+6 ' ) ||/(|u| 2 u) \\ X i.-»'({0,6]xM x )- ( 6 - n ) 

This follows from the retarded estimate on the Duhamel term in X s,h spaces (See for instance Lemma 
3.2 of [20], section 2.6 of [38], or Proposition 2.11 of [9] in the context of compact manifolds). As a 
result, we will be done once we show that ||7 (|w| 2 w) ||xi,-f>' ^ ll-^llx 1 . 6 - This will in turn follow from: 

3 

\\I(uiu5u3)\\ x i,- b > < Y\_\\lui\\ x i,b- (6-12) 

i=l 

We will deduce (|6.12l) from the cubic estimate (|3.2ip which we recall for convenience: 

3 
i=l 

for all s' > 1/2. 

To prove (|6.12p from (|6.13p . we split into two cases: 

Case 1: First, suppose that P<^ 3 NUi — ut for all i = 1,2, 3, where P^. 3 NU t — P^/-~£<^ 3 N u i- As a result: 

3 3 

\\I(uiu5u 3 )\\ xl> - b > ^ \\u 1 u5u 3 \\ xl} - b , < Y\_\\ui\\ x i,t < \\Iui\\ x i,b 

i=i i=i 

where in the first inequality we used that m is bounded by 1, in the second we used the trilinear 
estimate (|6.13l) , and in the third inequality we used that m(k) ~ 1 when \J— A < N. 

Case 2: Now suppose that one of the Ui, say u± for definiteness, satisfies P<^2NUi — 0, then using 
equation (I6.7[) we have: 

3 

\\I(uiU2~U 3 )\\ xl ,- b/ < N 1 ~ 8 \\u 1 U2~U 3 \\ x , _b/ < iV 1 ~ s ||ui|| Xs , 6 ||M2|| X »,b||w 3 || Xs ,b < ]^[ | \IUi | | x i,b 

i=l 

where we used (|6.7I) for the first inequality, the trilinear estimate (|6.13|) for the second inequality, and 
for the third we used (|6.7I) and the fact that ||ui||x 3i> ^ N ~ < " 1 ~ s * > \\I u i\\x 1 . b ■ 

Now (|6.12p follows by decomposing each m = P^ 3 NUi + (u, — P^ 3 NUi)- D 

Remark. The proof of the above proposition explains why we choose to work on the manifold M\ 
rather than M. Rescaling M to M\ allows for the normalization H/uqIIh 1 ~ 1 (cf. section IBTB"]) 
which yields for a time of existence 8 ~ 1. Without this rescaling, we would need to know the sharp 
dependence of 5 on ||7u||#i. Unfortunately, this sharp dependence is not provided directljQ by (|6.1ip 
and (|6.12p which lead to sub-optimal results. 



6.3. Decay of the modified energy. We are now ready to prove the almost conservation of the 
modified energy: 

Proposition 6.4. Lets > §, uq € H S {M\) withO < A < N. Ifu(t) solves (|6.2D and \\Iu\\ xl ,i+,r Q „ xM 
1, then 



though one can recover it by the scaling argument we do. 
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E[u(t)} - E[u(0)} 



< 



XNi 



3.14) 



for allO^t^S 



Proof. Computing the time derivative of E[u(t)] we get: 



d t E[u(t)] =d t E{Iu(t)) = Re / Iu t (~AIu + \Iu\ 2 Iu) dx 



Re / Iu t (-ilu t - AIu + \Iu\ 2 Iu) dx = Re / Iu t (~I(\u\ 2 u) + \Iu\ 2 Iu) dx 



Mx 



where we have used (|6.9p . Writing u(t, x) = ^2 n ^ n u(t,x) and using the fundamental theorem of 
calculus we get: 



E[u(t)]-E[u(0)] 

= Re>^ / / rn(ni)iT ni ut [771(712)771(71.3)777.(714) — /] n n2 u ir n3 u ir ni u dx dt. 
„. Jo JMx 



Since / is self-adjoint we can write this as: 



E[u(t)]-E[u(0)} 



Re 



Jm x 



m(ni) 



777(77.2)777(773)777.(77.4) 



7T ni (Iu t ) 7T„ 2 (Iu) 7T n3 (Iu) 7T n4 (Iu) dx dt. 



Using (|6.9p once more we get that: 



E[u(t)] - E[u(0)] = Re i (Termi - Term 2 ) 



(6.15) 



where 



777(771) 



777(712)777(773)777(774) 



7T ni (A/7t)7T„ 2 (Iu)TT n3 (Iu)lT ni (Iu)dx dt (6.16) 



Term 2 = £ / / 

„ JO J M x 



1 - 



777(774) 



777(772)777(773)777(774) 



?T ni (I (\u\ 2 u))n n2 (Iu)Tr n3 (Iu)Tr n4 (Iu)dx dt (6.17) 



6.4.1. Bound on Termi. We start by estimating Termi. Our goal is to prove that: 



Terrm <||/„|| 



1 



(6.18) 



For this we break u into a dyadic sum u = un where N = 2 J , for j = 0, 1, 2, . . . with un = P[n,2N) u 
for N > 1 and ui = -P[o.2)U. This reduces estimating Termi into that of an integral involving dyadic 
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frequency pieces at scale Nj . We will be able to sum over all those dyadic pieces by making sure that 
our estimates include a geometric decay factor in the highest frequency involved, that would allow 
using Cauchy-Schwarz to recover ||u||x»,i>- We present the details: 

First notice that: 



|A7tt|| x - I( i/2+ «S ||/u|| x i,i/2+. 



Therefore, in order to estimate Termi and prove (|6.18|) . we only need to show that: 



E 



N JO JM X 



1 - 



m(ni) 



771(712)771(713)771(714) 



1 



1 (4>l ) Kn 2 {4>2 ) 7r n 3 (03 )l"n 4 {4>i ) dxdt 
4 

n 



(6.19) 



< ^753(^1^2^3^4)° 1 1 ||<>,||v' 1 - 



i=2 



whenever the <f>j are spectrally supported on m ~ A/q Noticing that the above is symmetric^ in 
712,713,714, we will assume without loss of generality that 



N 2 > N 3 > N 4 . 



(6.20) 



We now conduct a case by case analysis depending on how N2 compares to Ni and N. In the particular 
case when M = T 2 or S 2 , one can directly assume that Ni < N2 because the right hand side of (|6.19[) 
vanishes if tii ^712 + 713 + 714 thanks to the sharp spectral localization mentioned in section U on those 
domains. Fortunately, this is a minor issue and will be dealt with using the following lemma: 



Lemma 6.5. There exists C > such that, if N± CN2, then for every q > 0, 



1 4 

L.H.S. of §M < q J] ||^|| X a, 1 / 2+ 
1 .'— 1 



(6.21) 



Proof. The proof of this lemma is straightforward using lemma I4~T1 We include it here for completeness. 
We will use the convention that for any function / : M\ — >• C, we shall denote by / : M — > C the 
pull-back function defined for y 6 M as f(y) = /(Ay). Rescaling the integral in the L.H.S. of (|6.19[) 
back into an integral over M and using the fact that Tv ni f(x) — 7TAn 4 /(f) w e get: 



L.H.S. of (IBTTgi) =A 2 



77i(t1i) 



771(712)771(713)771(714)/ Jo J M 



<A 2 V 1 



ni~Ni 



77l(71i) 



771(712)771(713)771(714) 



t 4 

2 (M)dt 

i=l 



^^In one particular case (Case 4 below), we won't prove the exact form in J6- 196 but rather use Cauchy-Schwarz in 
-/V2 (since we will have that Ni ~ N2) and geometric series summation to sum in ./V3 and JV4. 

19 Strictly speaking the symmetry is broken due to the existence of complex conjugates. However, these will not affect 
the analysis in any way and hence the treatment of the other cases is similar. 
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by applying lemma 14.11 Using the fact that ^ 1 since m is non-increasing and the fact that 



/ "r < 1 for j = 3,4, we can crudely bound the whole multiplier 1 -, — ^1,1, — r < n?. Also, 

by applying the Weyl asymptotics: 



#{^ : N } ^ ^/v ^ 2Nj with v an eigenvalue of -A g } < Nf 



and hence by Cauchy-Schwarz we get that Yl n -~N- I I^Ani 0i (*) 1 1 i 2 (A/) ^5 (^^i)H0j|U 2 (M)- As a result, 
we get: 



L.H.S. of @39D < — ^J* H(m)\\M^M)dt< Ap _4 f -6 /nil^WH^ 2 ^)^ 

^ — II 01 (*) IU? >a! ([0,i]xM A ) II 02 C*) II Ll x ([0,t]xM x ) 11^3 (*)IUf£|([0,t]xAfA) 1104 (*)|| ££([(),{] xM x ) 



< 

~ ATP 



4 

i\\X0,l/2+. 



□ 

Since the multiplier on the LHS of (|6.19l) vanishes when A2 <C Ai ^5 N, the above lemma is more than 
enough to give the needed decay on the RHS of (|6.19D . 

Remark. The same argument (with the roles of AT X and N 2 interchanged) gives that: 



1 4 

l.h.s. of < q w 11^11x0,1/2+ 



NX 
2 i=i 



unless N2 ~ maxjiVi, A3}. 

As a result, we will assume from now on that 

Ni < N 2 . (6.22) 

The analysis will be divided into several cases by comparing N 2 to N and the other frequencies: 

Case 1. N 2 <C N: The bound is trivially true since by (|6.22[) and (|6.20p the multiplier on the left 
hand side of (I6.19[) is zero. 

Case 2. N 2 > N A3 N4: This is the most delicate case. Applying lemma 15751 fand interchanging 
the roles of N± and A2) we may assume without any loss of generality that N 2 ~ N%. 

Split the dyadic interval [Ni,2N±) into J intervals I a of length A3 each, and the interval ^2,2^2) 
into K intervals Ip of length A3 each. Then J ~ K ~ As a result of this, we have: 
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L.H.S. of (I6J9I) 



E E 



JM X 



1 - 



m(ni) 



m(n2)m(n3)m(jii) 



Tni (0l)7Tn 2 (</>2)7Tn 3 (^3)^4 (<j>i)dxdt 



(6.23) 



We analyze (|6.23l) differently according to the respective locations of the two intervals I a and 1 p. To 
be more precise, we will call Si the sum in (|6 . 23[) when the interval I a is at distance at least 8-/V3 to 
the right of Ip (i.e. elements in I a are at least 8./V3 larger than those in Ip), S2 that when I a is at 
distance at least 8-/V3 to the left of Ip, and S3 when the two intervals are distances ^ 8-ZV3 apart. 



We start by estimating Si: First, we fix some notation. In this case N% ^ N%, hence we can write 
Ni = N2 + RN3 with R £ N U (since Ni, N2 are dyadic integers bigger than N3). As a result, 
[Nx,2N!) = lK _1 I q with I a = [N 2 + aN 3 ,N 2 + (a + 1)N 3 ). Similarly, [N 2 ,2N 2 ) = ifi^Ip = 



[N 2 + /3N3,N 2 + (/? + 1)^3). Since we are summing with (a, (3) £ Si, we have a - j3 > 8. 



E 

Si 



J M. 



= E^ 2 

Si 



m(ni) 



J m L 



m(ji2)Ta(n3)m{n\) 
m(ni) 



m(n2)°ra(n3)ra(nx) 



TTni (^O^ns (^2)^713 (03)7Tn 4 (4>i)dxdt 



TTni 01 7T„ 2 </>2 7Tn 3 </>3 ^4 4>±dydt 



where, as before, for any function / : M\ — > C, / : M — > C is given by f(y) = f(Xy) for every y £ M. 
Using the fact that ir n f(x) = ir\ n f(x) from (|2.1[) . we get by applying Theorem 14.21 that the latter 
expression is equal to: 



s JO JM X 



m(ni) 



E 

Si 



1 



771^2)111(713)171(714) 
77i(n\) 



771(712)771(713)777(714) 



TTni (4>l)^n 2 (02)7Tn 3 (03 ) Tti 4 (<f><l)dxdt 

X 2 (-2) 1 



((Atii) 2 - (An 2 ) 2 - (An 3 ) 2 - (An 4 ) 2 )' Vo -/m 



I'm 01 



x -B(7r n2 2 ,7r n3 03,7r„ 4 04) + C(7r„ 2 02,7rn 3 03,7rn 4 04) 



where 

ft) = O t+J+k= 2i (V7 * V J 5 * V k h) (6.24) 

<?«(/, 3, h) = O i+j+k ^(i-i) (V a i? * V7 * V J '<? * V*ft) . (6.25) 

In particular, both i? and C are trilinear operators that are linear combinations of products of differen- 
tial operators on M of the form Qi(f) Q2(g) Q3Q1) whose respective orders i,j, k satisfy i + j + k ^ 21 
and ^ k ^ I (where I will be chosen large enough so that the decay factor in \a — j3\ coming 
from the denominator n\ — n\ — n\ — ng would cancel a growth factor in \a — f3\ that comes from the 
multiplier thus giving a summable contribution (cf. (|6.30[) )). Redoing the scaling to go back to M\ we 
get: 
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E 



JM 



1 - 



m(ni) 



Si 



m(n\) 



771(712)771(713)771(714) 



777(772)777(773)771(774) 
(-2)' 

(nf - n2 - 77§ - nl)' ./ 



TTni (0l)7Tn 2 {^2)^3 OAs)?^ {4>i)dxdt 



M A 



Trsi {4>1 ) [Q2 (7Tn 2 02)<?3 (7Tn 3 03 ) 



where I.e. is a shorthand for the statement "a linear combination of and Qi(f)(x) 
differential operators of the type discussed in Corollary [ 

Let us denote: 



}i{TT ni cj>i)\dxdt. 

(6.26) 

:= [Q(/)](§) are 



777(771. n 2 , 7i 3 , 774) 



m(rii) 



(-2)' 



(77? _ n \ - n | _ n |y ■ 



777 (?72 ) 777 (773 ) 777 (774 ) 

Then the sum of Si can be written as a finite 0(1) linear combination of: 



(6.27) 



V] 777.(77,1 , 77 2 , 77 3 , T7, 4 ) / / n ni ((p 1 )Q2(irn 2 4>2)Q3(Trn 3 4>3)Q4('Kn 4 4 l 4)dxdt. (6.28) 
^ s JO J M A 

The estimate for this sum is obtained in two steps: 
Step 1: 

The first step is to apply a similar analysis to that in section [S] to bound the above integral using 
estimates on the multiplier 777. In fact, we will be able to write (|6.28[) as follows: 



T2I E E M^---,0a) I I Pi a <f> e 1 1 Q2(Pi^ e 2 2 )Q 3 (<t> e 3 3 )QM i )dxdt (6.29) 

A I,, fl ^o a JO JM X 



(a,/3)eSi 0i£ 

where A{9\, . . . , #4) is a summable sequence with 



J2 \A(e u ...,e i )\<(a-p)^- 1 



(6.30) 



and is a frequency modulation of (f)j (see (|6.35[1 ). In particular, (j)® 3 has the same X 5 ' 6 norms and 
frequency support as <f>j. 

Step 2 In the second step, we prove that for each fixed (#1, . . . , 84) € (Z/4) 4 , we have: 



1 

A 27 



JM X 



Pi a 4>? Q2 (02 2 )Qa W3 3 )^4(04 4 )d* * 



x ll^ Q 0l||x°, 1 /2+||-Pf /3 02||x o ' 1 /2+ 1 1 <?^3 1 1 J>CO. - 



(6.31) 



i=3 
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Combining ([5^5]) , (|OD|) , and (|QT) . one gets: 



(\ (N N "I 1 / 2 4 

^ 14(01,... A)l (A^ r |]Pj>lllxO^ + l]fi^2l|x°.V2 + nil^lixO.V 2+ 

<v-, al m (n 2 n 3 ) 1 (iy 3 iV4) 1/2 „ p , „ llp , .. A IU „ 

$ >. «- P TT 77; ^ ; ^ \ ll p r Q < PillAo.i/2+||n 3 2 ||x°-i/2+ I <P» x".V2+ 

^2 (N 2 N 3 (a - fi)) L A 11 

4 



: -"^11^9111x0. 1 /2+||n /3 02||x°. 1 /2+ I 1 11^411^0.1/2+ . 



By taking / ^ 3, we notice that by Schur's test (for example) and the fact that a (3 + 8 that 

™ < ^ (iV3 ^ 4)1/2 ]Ili^lix°.v 2+ (6.32) 

z— 1 

which will be enough to prove (16.191) . We now turn to proving Step 1 and Step 2: 

Proof of Step 1: The analysis here is almost the same as that in the proof of lemma IQl except for 
the localization to the intervals I a and Ip and the presence of the differential operators Qi . We include 
the details for the convenience of the reader. 

As in the proof of lemma IQ since m G I a = [N 2 + aN 3 ,N 2 + (a + 1)N 3 ), n 2 e Ip = [N 2 + (3N 3l N 2 + 
(,6 + l)N 3 ), n 3 ~ N 3 and n 4 ~ jV 4 we can write0: 



n\ = N 2 + aN 3 + N 3 ni with r?i € [0, 1) 
n 2 = N 2 + BN 3 + N 3 n 2 with n 2 G [0. ll 
n 3 = N 3 (l + n 3 ) with n 3 G [0, 1] 
n 4 = N 4 (l + n 3 ) with n A G [0, 1]. 

We now define the function \I/ : [0, l] 4 — > K given by: 

n 2 , n 3 , ra 4 ) = m(ni, n 2 , n 3 , n 4 ) (6.34) 

with rij given in (|6.33l) . Extend $ to a C°° compactly supported function on [—2, 2] 4 and then as a 
4— periodic function on R 4 . This allows us to express it in Fourier series: 

*(ni,...,fi4)= J2 e j(ei " 1+ - +94 " 4) A(6» 1 ,...,6i 4 ) 

6>iGZ/4 

with 



J2\ A (0i,---,6^\<\\ncm Ol mE 



2 ^with the obvious modifications if N3 or N4 is equal to 1, in which case we write ra 4 = n 4 £ [0, 2). 
21 See for instance [21] Theorem 3.2.16. 
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Notice that this already gives (|6.29p . In fact, 



(|6.28p =--Lt- ^2 ^2 m(ni,n 2 ,n3,n 4 ) 



where 



n 3 ~Ar 3j n4~iV4 



7T ni 01 Q2 (7Tn 2 02 ) Q3 ("Vi 3 03 ) Qi {^m 04 ) 



J Mi 



i E EE ^-a) 



A 2i 



(a,/3)GSi nie/ a ,n 2 e/ / 3 0;6Z/4 



7r ni 0i Q 2 (e iflafia tt„ 2 02 )Q 3 (e»3fts ^„ 3 4,3) Q 4 (e**™* 7r„ 4 4 )dxdi 



JM A 



4r E Y, M8i,.-.,9i) I I PiA 1 Q2(Pie4 2 )Q3(^)Q4(4 i ))dxdt 



E 

711 

E 



e ' ' 7T ni 01 



ifl 2 (r» 2 -(iV 2 + /3iV 3 )) 

e "3 7r n ,02 



E -ifl 3 ^3-JV 3 ) 
e x 7r n3 3 



2 y e ™ 4 7r™ 4 04- 



(6.35) 



In order to prove (|6 . 30[) and finish the proof of Step 1 all we need to do is prove that: 



H*||c=([o,i] 4 ) ^ 



N 3 



'N 2 (N 2 N 3 (a - (3)) r 
This estimate follows by direct verification using the following facts: 
(i) By the mean value theorem we have that 



m(ni) 



m(n 2 ) 



(ii) Similarly 



IP 



d(ni,h 2 y 



( _ m(N 2 + aN 3 + N 3 fn) \ 
V 771(^2+^3 + ^2)7 



<(^)H«( a -y3) 



N 2 



for any multi- index 7 with I7I ^ 2. This follows from the easily verified fact that: 

d l m(0 K 1 



MO ~ lei' 



for any (eM and / e N (see fl?TBJ)) 
(iii) A calculation shows that: 



(6.36) 



(6.37) 



(6.38) 
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n\ - {n\ + n\ + n\) =2N 2 N 3 (a - p + (n x - n 2 )) + (a 2 - (3 2 ) N 2 + iV 3 2 (n 2 - n%) 
+ 2N 2 (ani - (3n 2 ) - N 2 (l + n 3 ) 2 - iV|(l + h 4 ) 2 
= (2N 2 N 3 {a -13) + (a 2 - /3 2 ) N 2 ) Mni,n 2 , n 3 , n 4 ) 

with ip 2 : [0, l] 4 — ► R being C°°, bounded below by 1/2 (since in the Si sum we are estimating 
a — (3 ^ 8), and bounded above along with its derivatives by O(l). 
This gives that 



d{hi, n 2 , n 3 , n 4 p \ (n\ — n| — n| - n|)' 



< 



1 



(27V 2 iV3(a - 0) + (a 2 - /3 2 ) N 2 ) 1 
1 



{N 2 N 3 (a - /?))' 

Combining the above three facts one gets (|6.36l) . 

Proof of Step 2 The proof of (|6.31[) is an easy consequence of Corollary \3J>\ In fact, 



l.h.s.([OU ^ll^'Q^ll^o,^ 

(N^N^ 1 / 2 4 

T x \\Pl a ( f ) l\\x a - 1 / 2 +\\Plfi ( l ) 2\\x - 1 / 2 + Y\. II^Ha°' 1 /2+ 



i=3 



where we used (|3.18|) in the second inequality along with the fact that deg(Q2) + deg(Q3) + deg(Q4) ^ 21 
and ^ deg(Q 2 ) < I for i = 2, 3, 4. 

Similarly, one gets the same bound for S2 just by interchanging the roles of N\ and N 2 above. The 
bound for S3 = {I a ,Ip are at a distance ^ 8-/V3 apart} (near diagonal terms) is simpler. For each Ip, 
let S3 be the set of I a intervals that are at a distance less than ^ 8A^ from I p. Clearly there are ^ 19 
elements in S- 



0. 



E 



S3 



JM 



1 - 



m(n\) 



m(n 2 )m(ji3)m(n4 : ) 



EE E 



JM X L 



1 



m(ni) 



m(n 2 )m(n 3 )rn(n4) 



The sum 



E 



1 



m(ni) 



JM X L m(n 2 )m(n3)m(n 4 ) 
is bounded using lemma UTTl In fact, with 



7T ni (</>l)7T„ 2 (4>2)^n 3 (</ , 3)l"n 4 ((j>i)dx dt 



iV2 
iV 3 



1 



m(ni) 



771(712)771(77,3)771(77,4) 
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condition (I5.3|) was already verified in (|6.37|) and (|6.38p . The second condition (|5.4[) in lemma 15.11 
follows from (I3.3[) since: 



/ fl(t,x)f 2 (t,x)f 3 (t,x)f A (t,x)dxdt < ||/l/4|| i 2([0 ) i] x M A )ll/2/3|U| 3; ([0,t] X M x ) 
JMx 



(iV 3 7V 4 ) 1/2 



4 



— ^ ]J H/i|lxO.V2+ 

i— 1 



whenever /i, . . • , /j have the same frequency localizations as those under consideration. Therefore, 



S 3 JO ^Af A 



m(ni) 



m(n2)m(n3)m(n4) 



TTni (0l)7Tra 2 ((foVna (03 ^ru {<j>i)dxdt 



7V 3 (iV3iV4)V2 * 

$ 2^ at \ ll p f«'? , illjff.i/2+||n (3 02||xo-i/2+ _[_[ 



(6.39) 



i=3 



< 



N 3 {NzNj) 1 / 2 
N 2 A 



n 

i=l 



where we applied lemma 15.11 in the first inequality and used Cauchy-Schwarz in the Ip sum and the 
fact that | S3 I < 1 in the second. 

Combining (|6~32)) and (j6~39| we get: 



L.H.S. of « A II II^H — 

2= 2 

I 4 

~ Ajvi/2-jv° + ||0i||x ' i ' i/2+ n H^H* 1 ^ 

2 z=2 



which is what we want. 
Case 3. N 2 ~ N 3 > N 

The bound is obtained by invoking lemma I3~T1 In fact, by (|3.4p we have: 



/ fl{t,x) . . . f 4 (t,x)dxdt\ < \\f2h\\Ll([0,5]xM)\\f3U\\L^ x aO,S]xM) 
JM 



< 

A J V A 



i\\ X 0,l/2 + 



(6.40) 



(6.41) 



whenever fi are spectrally localized to frequencies ~ iVj, Ni < N2, and -/V4 ^ 7V3. This verifies the 
conditional estimate of lemma with B = (jVl ^ )1/2 and Y = X°^/ 2 +([0, 5} x M). 
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The symbol type conditions (|5.3[) are satisfied by the multiplier: 



m :- 



m{N 1 ) 



m(N 2 )m(N 3 )m(N i ) 



m(ni) 



m{n2)m{n'i)m(n4) 



since 



m(ni) 



777,(712)777(713)771(77,4) 



< 



m(ni) 



m(iVi) 



m(ri2)m(n3)m(n4,) m{N 2 )m{N 3 )m{N4) 



This follows from the fact that m y 11 ' > l since N 2 > Ni and m is non-increasing. Estimates of the 
derivatives follow as well since m itself satisfies the symbol-type estimate: 



dt l [, ~ tf 



which is due to (|6.6[F^ . As a result, 



3.42) 



L.H.S.of HOB < 



77l(-/Vi 



^^TTiuii 

"XX \\ < Pi\\x - 1 /"+([0,S\xM) 
i=l 



< 



m(N 2 )m(N 3 )m(N i ) A 

m(7Vi) JVi (7Vi7V 4 ) 1/2 

'771(^2)771(^3)771(^4) N 2 N 3 N 4 A 



^lllX-i.i/a+dO.^xM) H \\ ( Pi\\X 1 : 1 / 2 +([Q,S]xM)- 
i=2 

(6.43) 



We will estimate 



m(jVi) 

m(A r 2 )m(A r 3 )m(A r 4 ) 



depending on the relative position of iVi relative to TV: 



Case 3a N 2 ~ A3, iVi ^ A^: In this case, we use the explicit expressions for m(Ni),m(N 2 ),m(N 3 ) to 
get: 



1 ;yi— s iy— (i — s ) ,y 3 / 2 4 

L.H.s.of(Ei3<-— — — — ? f ir — -— 1 N>1/2 H0i|lx-w+nil^H^ a+ 

I iv 1/2+s 4 

2 v 47 4 „ (6.44) 

II 



<_ 



AjVl-aJV 2 a - 1/2 /2+([0 ' ilXAf) 

1 



II^IIa' 1 ' 1 /2+([o,5]xM) 



i=2 

4 



-^1/2-^0+ H^lllx-i.i/2+([0,5]xM) II H^Hxi.i/2+([0,5]xM) 
2 2=2 

where we have used the fact that N 2 ~ and JVi < in the second and third inequalities respectively, 
and that m{Ni)N\ > 1 for s > \ in the third inequality as well. 

Case 3b: N 2 ^ N 3 > N , Ni < N In this case, bound m(Ni) by 1 and use the explicit expressions for 
m(N 2 ) and m(N 3 ) to get: 

Here the key estimate goes as: 



22 As mentioned in the footnote to lemma [5.1l one only needs to verify the symbol estimates for two derivatives only. 
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L.H.S.of HPS! <- 



N: 



3/2 



A AT2(i- S )iV 2 " (1 " s) iV 3 " (1 " s) N 2 N 3 m{N 4 )Nl /2 



i=2 

4 



< 



< 



< 



N 



3/2 



1 iV 2 



(6.45) 



AJV 1 / 2 iV 57 I^HIX-i.i/a+dO^xM) nil^'Hx M/2+ ([0>]>!M) 



i=2 
4 



AA rl /2-iv 0+ l|(?;i|lx_1 ' 1/2+ ([ . <5 ] xM ) II llWIIxi.va+ao.^xM) 

2 z=2 



where we have used the fact that JVja ~ N3 and Ni < N in the second and third inequalities respectively, 
and that rn{Ni)N\'' 1 > 1 for s > i in the third inequality as well. 

Case 4 N 2 ^ N 3 > N: 

Recall that in this case Ni ~ iV-j. We will obtain the bound by invoking lemma 15.21 and the same 
bound for the multiplier as in Case 3. However, in this case, since N3 <£L Ni, N3 should replace Ni 
on the left hand side of f)6.41[) : 



L.H.S.of (jfU9)l 



<- 



miNx) (JV 3 iV 4 ) 1/2 A 1U „ 

■ 11 \\9i\\x . 1 / 2 +([0,S]xM) 



'm(iV 2 )TO(7V 3 )m(iV4) A 



1 



A m(N 3 )N^ /2 m(N 4 )N^ 

1 1 

'A ' 



< 



T7 2 "ll < / , l|lx- 1 - 1 /2+([o : 5]xAf) Y\. ll^llx 1 - 1 /2+ 



([0,5] xM) 



\-J^\\4 > l\\x^ 1 - 1 / 2 +([0J]xM) Y\_ \\ < t ) i\\x 1 . i / 2 +([0,S]xM) 



i=2 



m{N 3 )N 3 

~T nr" rll0illx- i > i /2+([o, ( 5]xM) Yl 

X N 1 - S N Q 2 ( s) • - 



(6.46) 



I >i\\X 1 - 1 / 2 +({0,5}xM) 



i=1 



< 



X jyi /2- iv°+ 1 1 91 1 1 * -1,1/2+ (IM x M ) II llWIIX^Va+do^xAf) 



where we used the fact that ~ Ni in the second inequality, that m{N4)N 1 / 2 ^ > 1 (since s > 1/2) 
in the third inequality, the explicit formula for m(N3) in the forth, and finally the fact that N 3 > N 
in the fifth. 

In this case, we do not have an exponential decay factor of the form iV 2 + in the denominator. However, 
since N% ~ N 2 on can use Cauchy-Schwarz to sum in Ni and N 2 and then use the geometric decay 
factor iV°+ to sum in N 3 and -ZV4. 

This finishes the proof of (|6.19p and hence the bound on Termi 



6.5.1. Bound on Terni2.' We now turn to proving the decay estimate for Term 2 in (|6.17[) . We will be 
able to prove better decay for Term 2 than that for Termi . In fact we will show: 

Term 2 < -J—\\I u \\ 6 xl , 1/2+ . (6.47) 
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This decay is proved by exploiting the high multiplicity of factors more than the exact shape of the 
multiplier (though of course the existence of the multiplier is crucial to restrict attention to the case 
where at least one factor has high frequency). In fact, since the specific form of the multiplier is not 
as important in this case, we will rewrite Term.2 in a form that is more convenient to do a multilinear 
analysis in LP spaces: 



Term 2 =]T f f 

„. JO JMx 



1 - 



ro(rii) 



771(712)771(713)771(714) 



JM x 



I (\u\ 2 u) Iululudxdt 



JMx 



7T ni (I (\u\ 2 u))ir n2 (Iu)n n3 (iu)7T„ 4 (Iu) dx dt 
I 2 (|7t| 2 u) uuudxdt. 



We now write u — un in Littlewood-Paley pieces as before to get: 



Tcrm 2 = E ( / / 1 (v>Ni 

N . e2 H \Jo JMx 



un 2 u n 3 ) Iun 4 Iun 5 Iu n 6 dx dt 



Jm 



I 2 (un 1 un 2 u n 3 ) un 4 un 5 un 6 dx dt 



(6.48) 



Let us denote by N max the largest of all Ni . . . Ng and N me d the second largest. As a result of 
the eigenvalue localization lemma |4"7TI and the fact that I 2 UNi = IuN t = ua^ if Ni ^ N/8, the 
contribution of the terms for which N max <C N to (|6.48|) is bounded by N~ q Yl i=1 WluWx 1 . 1 /'^ for 
some large enough q. In fact, writing unjU]^Un 3 = P^n {un!Un 3 un 3 ) + P>n (un/un^un-i), we get 
since I (P^nUn 1 u n 3 v> n 3 ) = I 2 (P^nUn^NsUns) that the contribution of the range N rnax <C N to 
(|jT4"5)l consists of: 



E 



IP>n {un 1 un 2 u n 3 ) un 4 un 5 un 6 dxdt 



I 2 P^ 



10 Jm x 



n (u Nl u N2 u N3 ) u Ni u Nro u Nf , dx dt 



Estimating \\IP >N (u Nl u N2 u N3 ) \\ L 2j M) and \\I 2 P >N (u Nl u N2 u N3 ) \\ L 2j M) by N c l\ l=1 \\u Ni \\ l%(m) 
for some large C and using Holder, a crude Sobolev embedding for the other terms, and the embedding 
_X*o,i/2+ c L^°L 2 one gets that the contribution of the range N rnax <C N is harmless. 

A similar argument using lemma 14.11 shows that the contribution of N max N me d is also harmless so 
we restrict attention to the case when N max > N and N max ~ N me d- 

Using the fact that / is bounded as an IP multiplier (cf. corollary 4.3.2 of 33 ), one can estimate the 
contribution of N max > N using Holder's inequality as follows: 



4 



Termf"""^ < ^ \\uN max \\ L lJ\u Nmed \\ L ^ *[[ IKvJU?,, 

N max ,N mei >N, »=1 



N 1/4 N 1/4 4 N 5/8 

^1/4^1/4 \\ UN m.^\\xo.i/2+\\UN mcd \\x'>. 1 / 2 + [I TTJs I \ u N t I lx°. 1 /2+ 



< 

N max ,N mei >N 
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As a result, we have: 
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Term 2 <- ^ m rAr._..WrAT_. J J |/ujv *»" llxl - 1/a+l|/ujv ^ llxl ' 1/2+ 



A ^ . m(N max )m(N med ) 

iV max £2 N , JV med ,Ni JV max 



ni V ■ 
-77rrll-f«jv < lljfi.v>+. 
i=1 m(Ni) 

Since A ma:E and A me d are > N, it is easy to see that for a > 1 — s, N!^ ax m{N max ) > N a and 
N^ led m(N mei i) > A". Similarly, for any fc > 0, m(k)k^ > 1 if /3 > 1 — s. Applying those estimates 
with a = 3/4- > 1/3 > 1 - s and /3 = 3/8- > 1/3 > 1 - s, we get: 



Term2 max ~ N < XN l /2 _ N max\\ Iu N mM \\x^+\\ Iu N mBd \\x^+ II H /u JVi llxi.V2+ 



1 HAfll 6 < 1 || 7-,.||6 

~AA 3 / 2 " 1 11 1* 1 ' 1 '** AA 1 - 11 11x1,1 J 



since A < A. □ 

6.6. Polynomial bounds on E[u(t)) and global well-posedness: We are now ready for the final 
step of the argument. Suppose that Uq E H s (M), then consider the function u$ : M x — > C given by 
u Q (x) = jU (f). Then \\u \\h s(Mx) = j?\\U \\ (m) and 



E[uo\ =\ f \VIu \ 2 dv x + - [ \Iu a \ 4 dv 

1 J Ms 4 J Mi, 

Using the fact, that uq{x) = jUq(j) = j ^2 TT Uj Uo(f), we calculate: 



\VIu \ 2 dv x =(-A x Iu , Iuo) l*(m x ) = {-A x I 2 u a ,u Q ) L 2 {Mx) 

2 __ 



AM °IaVAV; H"5 uo "iW-Z.^" l » m J m^^oii L 2 (m) 

where we have used in the last inequality that 1 1 7r zio 1 1 ^2 ( j\,^ A ) = IK^-^oll/^^) which comes from 
(|2.1I) and the fact that tt^ x2 uo(x) = j^uUo{j)- Splitting the sum in two cases Vj ^ 2(AA) 2 and 
Vj ^ 2(AA) 2 get that: 

N 2(l-s) 

^ \\ U 0\\ H *(M)- 



\VIu Q \ 2 dv x < -— _||t/ " 2 



'Ma 

Using the Gagliardo-Nirenberg inequality (which is scale-invariant) we get: 
(the scaling leaves the L 2 norm dimensionless). As a result, one gets: 



23 Here we used the second part of 1 12.711 corresponding to the case when Ni ^ A (i = 3, . . . , 6). The cases when Ni A 
are treated similarly and yield the bound Al / 2 j^-3/2- ^ ajv 1 - smcc ^ ^ 
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E [ u o] < ^^ll^o|| 2 ffs( M)- (6-49) 

Therefore, choosing A ~ ensures that E[u(0)] ^ i. By Proposition 16.21 there exists a S > such 

that: 



We can keep on using Proposition 16.21 to iterate the above inequality as long as E[u(kS)] < 1, and 
hence we can repeat the above procedure ~ AiV 1 / 2- times. Given T > 1, we choose N — N(T) 3> 1 
so that 

T - —^ N ^-. (6.50) 

Since the exponent is positive for s > |, N 3> 1 is defined for all times T> 1. Moreover, for ^ t ^ T, 
we have: 



\\U(t)\\ H s [M) = X s \\u(\ 2 t)\\ H s {Mx) < X s E[u(X 2 t)}^ 2 < X s ~ N^'WUqWh^m) 
since for ^ t ^ T we have E[u(X 2 t)] < 1. Using (16.50^ we get: 

TOIIff.(AO ^ T^-+\\Uo\\ H . {M) (6.51) 
finishing the proof of Theorem 11.31 

□ 

7. Spectral localization: Part II 

In this section, we give a proof of the results of Section 0] Here M d is a compact Riemannian manifold 
(possibly with boundary) with Riemannian metric g a p and A = A g is the Laplace-Beltrami operator. 
In what follows, we assume that are eigenfunctions of Laplace-Beltrami operator with Dirichlct 
or Neumann boundary conditions corresponding the the eigenvalues nf. The calculation leading to 
the proof of Theorem 14.21 and Corollary 14.31 is based on Bochner-type calculations involving the Ricci 
commutation identities used to commute covariant derivatives. This is basically the way to generalize 
some integration by parts manipulations to the case when the involved functions are contracted tensors 
rather than just functions. We will do calculations without resorting to a prefered coordinate system. 
The notation we will use is fairly standard: We use abstract index notation^!, repeated indices are 
summed, and g is used to raise and lower indices. We will quickly review some elementary concepts 
from Riemannian geometry for the sole purpose of fixing notation. For a more comprehensive treatment 
we refer the reader to any of the standard texts on the subject 32] [TO] or the first chapter of [T7] for a 
review of the commutation identities we will use. 



24 in particular the indices do not correspond to any preferred coordinate system. 
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7.1. Review of some elementary Riemannian geometry concepts. We denote by V the Rie- 
mannian connection associated to g. This is the unique connection that is torsion free and for which 
the metric g a p is parallel. In other words, 

Vg Q/ 3 = (Metric Compatibility) 
and for any vector fields X and Y: 

V X Y - VyI = [X, Y] (Torsion Free) 
where [X, Y] = XY - YX is the Lie Bracket. 

Ricci commutation formulas: For any C°°(M) function /, it follows from torsion- freeness that: 

V Q V /3 / = V/jVcJ. 

Two covariant derivatives only commute when acting on functions. In particular, if one takes two 
covariant derivatives of higher order tensors (for example taking two covariant derivatives of V/), they 
need not commute. This leads to the definition of the Riemann curvature tensor R which acts on 
tensors to measure the error incurred from commuting two covariant derivatives. We start with the 
case of a vector field Z, in this case the Riemann curvature tensor is defined as: 

R{X,Y)Z := Vx(VyZ) - V Y {V X Z) - V [xy] Z 
or more simply using abstract index notation: 

iW^ 7 = (VaV^ - v v a )z s E 

What will be important to us is that R is actually a tensor, in the sense that Z is not differentiated 
when computing R(X, F)iQ 

The Riemann curvature tensor R s a a does not only give the commutation rules for covariant derivatives 
acting on vectors, but also for those acting on general (fc, I) tensors T — Tg^'"$ l . We first start with 
1-forms: when T — w 7 is a 1-form, since V ' g — 0, we have: 

(V a V^ - V^V«)w 7 = -Ri^us 
More generally we have for a (fc, I) tensor T = T^'"^ ■ 

(v Q v, - v^v Q )T£;;£< = £ R^ m T^r mii+ ^ 11 - E (7-1) 

i=l j=l 



The coordinate-free notation and the index notation are related by: 

(R(X, Y)Z) S = R aM 6 X a Y l3 Z'< 

for any vector fields X = X a , Y = , and Z = Z~* . 

2 ^This is in contrast to the fact that the commuting two pseudo-differential operators of order 1 gives a pseudo 
differential operator of order 1. 
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(see for example the first chapter of (T7J ) . As an application of the above Ricci commutator identities 
we have: 

Commutator of A and V: 

Recall that the (rough) Laplace operator (also known as the Bochner LaplaciarQ) is defined on general 
tensors as: 



A = divV = trace g V 2 = V Q V" = g af *V a Vp. 

Lemma 7.2. (i) For any function f 

(AV a — Vq.A)/ = Ric Q ^ V 9 / 

where Ric Q( g := R ia fp is the Ricci tensor. 
(ii) For any tensor T 

AVT - VAT = 0(R* VT) + O ((V Ric) * T) 

where, given any tensors A and B, A*B denotes some contraction of A®B and 0{A* B) 
denotes a linear combination of contractions of A®B. 

Proof. We present the simple proof as a warm up for the calculations to come. The proof of (i) follows 
easily from: 

AV Q / = VpV l3 V a f = V^V Q V' 3 / = V^V 3 / + %/V7 = V Q A/ + Ric Q/J V^/ 

where in the second equality we used that covariant derivatives commute on functions, while in the 
third and forth we used the definition of R and Ric respectively. 

The proof of (ii) follows similarly using (|7.1[) (see the first chapter of 1 7J for a more comprehensive 
review) . 

□ 

7.3. Proof of Theorem 14.21 We are now ready to present the calculation leading to (|4.3I) . For 

i = 1, ... ,4, let ej be an eigenfunction of the Laplace-Beltrami operator with Dirichlet or Neumann 
boundary conditions corresponding to eigenvalues — nf. Denote: 



e\(x) . . . e4,{x)dx. 

M 

Then by Green's theorem^ 

n\A = / (-Aei)e 2 e3e4(fa; = / ei(-A)[e 2 e3e4]da 

JM JM 

But 



This is in contrast with the closely related Hodge Laplacian which will not be of concern for us. 
"^Recall that Green's theorem states that f M (uAv — vAu)dx = fg^i 11 ^ ~ V T^0^ where d/dn denotes the normal 
derivative on the boundary and dS is the induced measure on dM. Since we are either assuming Dirichlet or Neumann 
boundary conditions, all boundary integrals vanish. 
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/ ei(-A)[e 2 e 3 e 4 ]<iz; = / e±(x) ((-Ae 2 )e 3 e 4 + e 2 (-Ae 3 )e 4 + e 2 e 3 (-Ae4)) dx 
hi Jm 



-2 e x (V Q e 2 V Q e 2 e 3 V"e 4 + e 2 V Q e 3 V"e 4 ) eix 

z ( n l + n l + nl)A Q -2 ex (V a e 2 V"e 3 e 4 + V Q e 2 e 3 V"e 4 + e 2 V Q e 3 V"e 4 ) dx. 
Jm 



As a result, we get: 



where 



A - 2Al 
^0 — 



nr — ns — n% — n. 



Ax = ex (V Q e 2 V Q e 3 e 4 + V Q e 2 e 3 V a e 4 + e 2 V Q e 3 V Q e 4 ) dx. 
Now we repeat the same procedure for Ax : 

71^!= / (-Aei)(V Q e 2 V Q e 3 e 4 + V a e 2 e 3 V Q e 4 + e 2 V Q e 3 V Q e 4 )cfe 
Jm 

= / ei(-A)[V«e 2 V a e 3 e 4 + V Q e 2 e 3 V Q e 4 + e 2 V Q e 3 V a e 4 ]da; 

JM 

by Green's formula. We now compute (— A)(V Q e 2 V a e 3 e 4 ): 

(-A)(V Q e 2 V Q e 3 e 4 ) = (-AV a e 2 )V Q e 3 e 4 + V Q e 2 (-AV Q e 3 )e 4 + V Q e 2 V Q e 3 (-Ae 4 ) 

- 2 (V Ql V Q0 e 2 V Ql V Q0 e 3 e 4 + V Ql V cto e 2 V Q °e 3 V Ql e 4 + V Qo e 2 V Ql V Q °e 3 V Ql e 4 ) 
= (-V Q A e2 )V Q e 3 e 4 + V Q e 2 (-V a Ae 3 )e 4 + V Q e 2 V Q e 3 (-Ae 4 ) 

- Ric af} V /5 e 2 V Q e 3 e 4 - Ric Q/3 V Q e 2 V /3 e 3 e 4 

- 2 (V Ql V Q0 e 2 V Ql V Q0 e 3 e 4 + V Ql V Qoe2 V Qo e 3 V Ql e 4 + V Qo e 2 V Ql V Q °e 3 V Ql e 4 ) 
={n\ + n§ + n\) (V Q e 2 V Q e 3 e 4 ) - 2 Ric Q(3 V' 3 e 2 V Q e 3 e 4 

- 2 (V Ql V Q0 e 2 V Ql V Q0 e 3 e 4 + V Ql V ao e 2 V Q °e 3 V Ql e 4 + V Qo e 2 V Ql V Q °e 3 V Ql e 4 ) 

(7.2) 

where we have used lemma 17.21 for the second inequality. Similarly, one gets: 



(-A)(V Q e 2 e 3 V Q e 4 ) ={n\ + n§ + n\) (V a e 2 e 3 V a e 4 ) - 2 Ric Q/3 V' 3 e 2 e 3 V Q e 4 

- 2 (V Q1 V QQ e 2 V Ql e 3 V Qo e4 + V Ql V Qo e 2 e 3 V Ql V Q °e 4 + V ao e 2 V Ql e 3 V Ql V Q °e 4 ) 

(7.3) 



and 
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(-A)(e 2 V a e 3 V Q e 4 ) ={n\ + n\ + n\) (e 2 V a e 3 V a e 4 ) - 2 Ric Q/3 e 2 V 3 e 3 V a e4 

- 2 (V Ql e 2 V Ql V ao e 3 V Q °e 4 + V Ql e 2 V ao e 3 V Ql V Q °e 4 + e 2 V ai V Qo e 3 V Ql V a °e 4 ) 

(7.4) 

Adding 0, CSJ), (HI), we get: 

n\A x = (nj + n\ + nl)A 1 - 2A 2 

where A 2 is of the following form: 



M ■= I ei (S 2 (e 2 ,e 3 ,e 4 ) + C 2 (e 2 , e 3 , e 4 )) da; 

JM 

where -B 2 (/, <7, h) and C 2 (/, g, h) are trilinear operators that can be expressed schematically as: 



B 2 (f, g, h) = O i+j+k = 4 (V7 * V J 5 * V k h) (7.5) 

C 2 (/, 5 ,/i)=Oi +i+fe<2 (i?*V7*V J .g*v7). (7.6) 
Now suppose, as an induction hypothesis, that: 

4 , — Z^L A 

s^n—l — 9 o 9 9 ^-n 

where 



A n = / e 4 (_B n (e 2 ,e 3 ,e 4 ) + C„(e 2 ,e 3 ,e 4 )) 



and 



B„(/, 5 , fc) = CWfc=2„ (V7 * V J <? * V7) (7.7) 

C n (f,g,h) = a i+J+fe ^ 2(n _ 1) (V a i? * V7 * V J 5 * V7) (7.8) 

where a is some exponent (which can be calculated explicitly in terms of n) signifying a number of 
derivatives applied to the curvature tensor. 

Then, as before, by Green's theorem: 



n l A n = I (-A)ei (£„(e 2 , e 3 , e 4 ) + C„(e 2 , e 3 , e 4 )) dx = I e 1 ((-A)B n (e 2 , e 3 , e 4 ) + (-A)C„(e 2 , e 3 , e 4 )) 
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and our goal is to write this as 



n\A n = (nl + n| + nf)A n - 2 / e x (B n+1 (e 2 ,e 3 ,e 4 ) + C n+1 (e 2 ,e 3 ,e 4 )) dx 

Jm 



with B n+ i, Cn+i as in (|7.7I) and (17.81) with n replaced by n + 1. For this we use the rules to commute 
A and V in lemma EHJii)- In fact, 



(-A)B n (e 2 , e 3 , e 4 ) = O i+i+fe=2 „(-A) (V*e 2 * V J e 3 * V fe e 4 ) 

= O l+3+k=2 n\ ((-A)V l e 2 * V J e 3 * V fc e 4 + V l e 2 * (-A)V J e 3 * V fc e 4 + V I e 2 V J e 3 (-A)V' c e 4 ) 

+2 (V a V ! e 2 * V Q V J e 3 * V fe e 4 + V Q V l e 2 * V J e 3 * V Q V fc e 4 + V'e 2 * V Q V J e 3 * V Q V fe e 4 ) 

- O l+J+k=2 n\ (V*(-Ae a ) * V J e 3 * V fc e 4 + V l e 2 * V J (-Ae 3 ) * V fc e 4 + V I e 2 V J e 3 V' c (-Ae 4 )) 

+2 (V a V i e 2 * V Q V J e 3 * V k e A + V Q V l e 2 * V J e 3 * V a V fc e 4 + V i e 2 * V^V^eg * V a V fc e 4 ) 

+O i+j+k<2n (W a R * V l e 2 * V J e 3 * V fe e 4 ) 



We should note that for the third equality above, we used lemma I7.2f ii) inductively applied to the 
tensor V i_1 e 2 (and similarly for V J ' _1 e 3 and V fc_1 e 4 ) to give AV(V i_1 e 2 ) - VA(V i_1 e 2 ) = 0(R * 
V J e 2 + V Ric *V l-1 e 2 ) which is of the form above. Doing the same thing for AV l_1 e 2 over and over 
we get the result claimed. As a result, we get: 



(-A)B n (e 2 , e 3 , e 4 ) =(n 2 + n 3 + n 4 )B n (e 2 , e 3 , e 4 ) 



+ O i+j+k=2n V Q V*e 2 * V Q V->e 3 * V fc e 4 + V Q V ,; e 2 * V J e 3 * V Q V fe e 4 

+ V l e 2 * V Q V J e 3 * V Q V fc e 4 ) + O l+]+k<2n (V Q i? * V l e 2 * V 3 e 3 * V fc e 4 ) (7.9) 
--(nl +nl + n\)B n (e 2 , e 3 , e 4 ) + O i+ j+fe=2( n +i) (V ,; e 2 * V J 'e 3 * V fc e 4 ) 

0^'ij',fe<n+l 

+ O i+j+k ^ 2n (V a i? * V*e 2 * V J e 3 * V fc e 4 ) . 



Obviously, the second term above will join B n+ i whereas the third will be part of C n +i- The compu- 
tation for AC,, is similar: 
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-A)C„(e 2 ,e 3 ,e 4 ) =e>; +J - +fe<2 (n-i) (-A) (V a i? * V 4 e 2 * V J e 3 * V fe e 4 ) 



=O i+i + fe <2(n-i) ( V a i? * (-A)V*e 2 * V J e 3 * V fc e 4 + V Q i? * V l e 2 * (-A)V J e 3 * V*e 4 
0<ij',fe<n-l 



V Q i? * V*e 2 * V J e 3 * (-A)V fe e 4 ) + O i+j+k ^ 2n (V a 'i? * V'e, * V J e 3 * V 



fe e 4 



=O i+i + fe <2(n-i) V a i? * V l (-Ae 2 ) * V J e 3 * V fc e 4 + V Q i? * V l e 2 * V J (-Ae 3 ) * V fc e 4 
+ V a R * V*e 2 * V J e 3 * V fe (-Ae 4 ) ) + O i+j+k ^ 2n (V a i? * V l e 2 * V J e 3 * V k e 4 ) 

by applying lemma \77I\ inductively as before. As a result, we get: 



(-A)C n (ea, e 3 , e 4 ) = (n| + n 3 + 7^)C„(e 2 , e 3 , e 4 ) + C^+^n (V a i? * V l e 2 * V J e 3 * V fe e 4 ) . (7.10) 



The last term will join C„+i(e 2 , e 3 , e 4 ) to give that: 



n x A n = (n 2 + n 3 + n 4 )A n + / e 4 (_B n+ i(e 2 , e 3 , e 4 ) + C n+ i(e 2 , e 3 , e 4 )) da; 

where B n+ i and C„+i are trilinear operators as in (|4.5I) and (I4.C5[) respectively, with n replaced by 
ri + l. This concludes the induction proof. D 

7.4. Proof of Corollary 14.31 : The proof of Corollary 14. 31 will follow by applying Theorem 14.21 and a 
variant of the bilinear eigenfunction cluster estimates of [S] [TU] after dealing with a couple of technical 
problems. Let us first recall the bilinear eigenfunction cluster estimates (also called bilinear Sogge 
estimates) from [TP] , 

Proposition 7.5. (Bilinear eigenfunction cluster estimates [10] J 

Let x G <S(K). For A G R, denote by x\ — x(V~A — A) i/ie spectral projector around A. -For any fi, 

IIxa/Xm5IU 2 (m) <i A(d,/x)||/|| I , 2 ( M -)||s'|| i2 (M) (7.11) 
for all /, g € L 2 (M) where A(d, /x) was defined in (|4.8[) . 

The proof in 10 is based on Sogge's parametrix representation of xa/ hi local coordinates (see |33|). 
More precisely, for every N ^ 1, one has the splitting: 

Xa/ = X^T x f + R x f 

with 

||-Ra/||h''(M) ^Ar : fe ^ _ H/IU^M) 

and in a system of local coordinates around each x$ € M , has the following parametrix representa- 
tion: 
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Txf(x)= [ e lX ^a(x,y,X)f(y)dy (7.12) 

where a(x, y, A) is a polynomial in A -1 with smooth coefficients supported in a compact subset {(x, y) £ 
V x V : V C K d (compact)} and —(f>(x,y) = d g (x,y) is the geodesic distance between x and y. By 
taking N large enough, the proof of (|7.11[) reduces to proving estimates for oscillatory integrals of the 
form CLU (cf. HD]]). 

As a result of this, one can directly deduce by the exact analysis as that leading to (|7.11[) that the 
following estimate holds: 

||(VW) * (V k x^)\\ L 2 iM) < \^ k A(d,ri\\f\\ LHM) \\g\\ LHM) (7.13) 

since any contraction of the form (V j Xa/) * (S k X^9) can be written using the above splitting of 
Xxf = ^~^~T\f + R\f and Xp.9 = + -^m as a nrLear combination of products of operators of 

the same form for which the analysis in [10] applies (what is crucial in this reduction is that R is a 
smoothing operator for large enough N and so it does not affect the analysis in any significant way. 
Also, a is just a symbol of order in a; and is a polynomial in j-, so any derivative falling on it does 
not affect the analysis in any way either). 

We now turn to the proof of proposition 14.31 Without loss of generality, we assume that / = 
1[a,a+i](V I: A)/, g = (\/ = A)fif, and h = l [i , )I/+1] (\/ = A)/i. Write 

/ h(x)f(x)g(x)dx= V] / n ni fnr n2 f ir n3 g dx 
Jm ni e[v, v +i] jM 

n 2 G[A,A+l] 

where n n is the projection onto the n 2 -eigenspace. By lemma l4~3l (with — l)we get: 



f (—2) f 

/ h(x)f(x)g(x)dx = ^ ~~, ; n m h { B (^n,f, n n3 g) + C(tt„ 2 /, ir n3 g) dx (7.14) 

Jm ni eW,u+i] {ni-nj-ni) Jm 

n 2 6[A,A+l] 
n 3 e[fi,fi+l] 

where 



B(n n J,TT n3 g) = (V j tt„ 2 /) * (V J 7r„ 3 /) 

and 



C(Tr n J,ir ri3 g) = J+K2( . 7 _ 1) (V°E* (V%„ 2 /) * (V j ir n3 g)) . 

The integral in (|7.14l) is similar to the ones studied in section [5] Essentially, one would like to argue as 
follows: since v = A + Kfi + 2, then n\ — n^ — n\ ^ 2K\(i and hence one would like to estimate (I7.14[) 
as follows: 
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LHS of (UP <„ 



(KXfj,)' 



h(x)(B(f,g)+C(f,g))dx 

M 



7 — (KXfJ,) J — ll ft ll L2 ( M )ll/ll£ 2 (M)||fl , ||£ 2 (Af) 



where the first inequality will be justified in what follows whereas the second inequality comes from 
(17. f 3p and the fact that all derivatives of R are bounded to get the second inequality. To justify the 
first inequality, we argue similar to what we did in section [5] We include the details for completeness. 
As before, we write ri\ = v + n, n-i = A + r 2 , and = /.i + r$ with (n, r 2 , r 3 ) € [0, l] 3 . As a result, 

2 J 2 J 

we get that 7-2 — t — jtt = tt — ; — 75 — ^ , ,, — 7 — : — c^tt- But 

3(xi, £2,273) : = r 

((^ + ^i) 2 -(A + .T 2 ) 2 -(M + a ;3 ) 2 ) J 

is a smooth function on [0, l] 3 that is bounded along with its derivatives by jj^jj^jj- In fact, this follows 
from the following estimates: 



2 J 2 J 1 
< = < 



(„2 _ n \ - n \y ^ ( v 2 _ (A + 1)2 _ ( M + !)2)J ((A + Kp + 2) 2 - (A + l) 2 + (/1 + f ) 2 ) J " (K\fi) J 

since if > 1. A similar estimate holds for the first two derivatives. Multiplying f(xi, X2, £3) by a 
compactly supported function on [—2, 2] 3 and extending the resulting function periodically to R 3 gives 
a 4— periodic function on R 3 . Expressing this function in Fourier series we get that E(a;i, X2, X3) = 
j: giem A(9 1 ,9 2 ,d 3 )e^+ e —+ 9 ^ with 

E 1^1,^^3)1 < (7-15) 

6>ieZ/4 v w 

since A{6\, 62, O3) are the Fourier coefficients of a C 2 function whose C 2 norm is bounded by rgx^p ■ 
With this in hand, we write: 



f h(x)f(x)g(x)dx = £ E MPufaMe 

M neS-i/n[o,i] e,ez/4 

r 2 es-An[o,i] 
r 3 es-/in[o,i] 



+02l"2 +#3»"3 ) 



x / ir^+^hiB^x+rJ^f.+^g) + C(ir\ +r2 f,'K u , +r3 g)dx 



M 



where we used 5 to denote the set {n e R : n 2 € spectrum of(—A 9 )} and 5 — a = {n 
n e S}. Letting h 9l = Er 1 es-,n[o ! i] e ' 9 "' 1 Vn'», fe 2 = Er 2 es-An[o,i] e ' if,2r27r A+r 2 /, and 
Er 3 es-/jn[o 1] el03T37r n+r39> we get, using the fact that -B and C are multilinear, that: 



/ h(x)f(x)g(x)dx = E ^1^2,^3) / h 8l [B(fg 2 ,g 03 ) + C(fe 2 ,gd 3 ) dx. 



But for each fixed (61,62, 63) € Z/4, we have the estimate: 



GWP OF CUBIC NLS 



47 



I he 1 (B(fe 2 ,ge 3 ) + C(fe 2 ,ge^j dx ^||VHl 2 (m)II (B(fg 2 ,gg 3 ) + C(fg 2 ,g e ^j \\l'(M) 

<X J fi J X(d, fj) \\h 9l | | L 2 (M) \\fg 2 1 | L 2 (M) | \gg 3 \\ L 2 [M) 
=A J ^ J A(d,/z)||/i|| L 2 (M) ||/|| i 2 (Jl/) || 3 || L 2 (M) 



where we used Cauchy-Schwarz for the first inequality, (I7.13|) in the second, and the fact that frequency 
modulation leave that L 2 norm invariant in the third. 

As a result we get: 



h(x)f(x)g(x)dx 



M 



<\ J ft J X(d, ft) Y, \ A ^1,02,O 3 

\e 4 ez/4 

\ J H J A(d, n) 



< 



HMU 2 (M)II/I|l 2 (M)II#IIz, 2 (M) 



L 2 (Af)ll/IU 2 (M)llfflU 2 (M) 



K j 



as desired. 



L 2 (M)\\j\\L 2 (M)\\g\\L 2 (M) 
□ 
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